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Abstract 


The study of the propagation of soliton waves in optical waveguides constitutes a fundamental 
and very interesting subject of research, particularly in the field of high-debit optical 
transmissions. In this context, this work focuses on the theoretical study of the nonlinear 
propagation dynamics of Bragg solitons in a fiber Bragg grating in femtosecond régime, 
where the propagation of the puises is modeled by a System of coupled mode équations 
(NLCM), more complicated with high order terms of cubic, quantic and septic nonlinearities. 
In order to solve this non intégrable System of équations, we hâve used the multiple scales 
analysis methods to reduce the System of coupled nonlinear équations (NLCM) into a well- 
known perturbed nonlinear Schrôdinger équation (PNLS) that describes the nonlinear puise 
propagation at the edges of the PBG structure. Then we applied the fonction of the combined 
complex amplitude, to solve this complicated équation which contains higher order nonlinear 
effects such as the quintic and septic non-Kerr terms, and higher order linear effects such as 
third- and fourth-order dispersions. The solutions obtained are Bright and dark Bragg solitary 
wave type. In addition, parametric conditions for the existence of these solutions hâve been 
presented. We hâve also calculated the important and interesting physical parameters such as 
power and puise width for Bragg solitons. 
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Résumé 


L’étude de la propagation des ondes de types soldons dans les guides d’ondes optiques 
constitue un sujet de recherche fondamental et très intéressant, en particulier dans le domaine 
des transmissions optiques à haut débit. Dans ce contexte, ce travail s’intéresse à l’étude 
théorique de la dynamique de propagation non-linéaire des solitons Bragg dans un réseau de 
Bragg fibré en régime femto-seconde, où la propagation des impulsions est modélisée par un 
système d’équations des modes couplés (NLCM), plus compliqué avec des termes d’ordre 
élevé de non linéarité de troisième, cinquième et septième ordre. Afin de résoudre ce système 
d’équations non intégrable, nous avons utilisé en premier la méthode perturbative multi- 
échelle qui a permis la réduction de ce système d’équations non linéaire couplées (NLCME) 
en une équation dite « équation de Schrôdinger non linéaire perturbée (PNLS) ». Cette 
dernière décrit la propagation non linéaire des impulsions aux bords de la zone du gap de la 
bande photonique. La résolution de cette équation compliquée contenant des termes non-Kerr 
et des termes de dispersion de troisième et quatrième ordre a nécessité l’application de la 
fonction de l’amplitude complexe combinée. Les solutions obtenues sont de type solitons 
noirs et brillants. En outre, nous avons également présenté les conditions d’existence de ces 
solutions ainsi que leurs importants paramètres physique tels que la puissance et la largeur de 
l’impulsion. 
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Introduction généra 1 e. 


La communication à distance a toujours été un souci majeur pour l’humanité et le 
développement des moyens et des techniques permettant le transfert d’informations n’a cessé 
d’évoluer dans le temps. Au début, l’homme ne pouvant utiliser que des signaux sonores et 
visuels a vu que ces techniques répondaient seulement à des besoins de voisinage et ne 
permettaient que la transmission de messages rudimentaires. Il a fallu attendre le 19 iem siècle 
avec l’invention du télégraphe électrique en 1837 par Samuel Morse, puis l’invention du 
téléphone en 1875 par Alexander Graham Bell qu’on a pu assister a la naissance des 
télécommunications modernes et depuis une vraie révolution a suivie en prenant de l’ampleur 
de jour en jour. 

L’utilisation des systèmes de transmission sur câble durant le 19 iem siècle a amélioré les 
performances du monde de la communication, mais ce qui a suscité un grand intérêt reste la 
recherche de la qualité et le débit de la transmission. Le développement des systèmes de 
communication à fibre optique au milieu du 20 iem siècle a complètement révolutionné ce 
domaine. En effet l’idée d’utiliser la lumière comme porteuse de l’information surgit de 
nouveau en 1960 avec l’invention du laser [1], Actuellement la fibre optique représente le 
meilleur moyen pour transporter de très hauts débits d’information numérique sur de longues 
distances cependant, elle reste confronté à deux types de limitations majeures qui sont 
l’atténuation et la dispersion. 

Ces limitations peuvent s’affranchir par l’introduction d’un nouveau concept « soliton », 
découvert pour la première fois par l’ingénieur écossais John Scott Russell en 1834 [2]. Le 
soliton optique a été découvert théoriquement par Hasegawa et Tappert en 1973 [3], puis son 
observation expérimentale en 1980 a été démontrée par Mollenauer et al [4]. Le soliton ou 
l’onde solitaire est un phénomène physique apparu dans les milieux non-linéaires et 
dispersifs, défini comme une perturbation localisée pouvant se propager sans déformation sur 
de longues distances, et ne perdant rien de ses caractéristiques ni par amortissement, ni par 
effet de choc éventuel avec une perturbation de même nature. C’est pour cela que la technique 
de transmission par soldons optiques est une alternative attrayante pour l’optimisation des 
performances et l’élimination des pertes de propagation. 
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Après la découverte du soliton optique, de nombreux laboratoires de recherche en 
transmission de l’infonnation ont travaillé sur la dynamique non-linéaire et la stabilité des 
solitons dans les guides d’ondes optiques tels que les réseaux de Bragg fibrés. Comme de 
nombreux guides d’ondes optiques, les réseaux de Bragg fibrés sont susceptibles de donner 
lieu à la propagation d’ondes solitaires. Ce phénomène a été mentionnée pour la première fois 
en 1987 par Chen et Mills sous la forme d’un soliton particulier appelé « soliton Bragg » [5]. 
Le soliton Bragg peut se propager sans déformation dans le réseau avec des fréquences située 
dans la bande interdite dite (Gap soliton), ou proches mais en dehors de celle-ci. 
Théoriquement, la propagation des impulsions optiques dans le réseau de Bragg est modélisée 
par les équations des modes couplées non-linéaires (NLCE) [6], et les solitons Bragg sont les 
solutions analytiques de ces équations. En recherchant des solutions situées exactement au 
centre de la bande interdite, Christodoulides et Joseph sont les premiers à avoir développé 
une méthode basée sur le modèle de Thirring massif [7]. Les travaux d’Aceves et Wabnitz 
[8] viennent compléter cette méthode et ont pu la généraliser pour n’importe quelle position 
dans la bande interdite. 

Récemment plusieurs études théoriques et numériques ont été réalisées pour clarifier et 
classifier les conditions de génération des solitons Bragg, et la dynamique non-linéaire de leur 
propagation. Pour cet objectif notre contribution a une vision purement théorique. Elle est 
axée sur l’étude analytique de l’existence des ondes solitaires dans un réseau de Bragg fibré 
en régime femto-secondes, où la propagation des impulsions est modélisée par un système 
d’équations des modes couplés comportant des termes d’ordre élevé de non linéarités de 
troisième, cinquième et septième ordre. Les applications physiques de ce modèle se situent 
dans un domaine de pointe qui est celui de la génération des solitons Bragg dans un milieu 
fortement non linéaire comme le verre de chalcogénure. Ce type de verre présente non 
seulement des non-linéarités de troisième et cinquième ordre mais aussi de septième ordre. 
Afin de résoudre ce système d’équations non intégrable, nous avons utilisé en premier la 
méthode perturbative multi-échelle qui a permis la réduction de ce système d’équations non 
linéaire couplées (NLCME) en une équation dite l’équation de Schrôdinger non linéaire 
perturbée (PNLS) [9]. Ensuite on a appliqué la fonction de l’amplitude complexe combinée 
[10] pour trouver les solutions exactes de cette équation. 
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Ce manuscrit est organisé comme suit : 

Le premier chapitre présente un rappel général sur la théorie de propagation des ondes 
optiques dans la fibre optique comme guide d’onde dispersif et non-linéaire. Cette 
propagation est modélisée par l’équation de Schrôdinger non-linéaire, qui peut facilement être 
résolue numériquement par la méthode Split-Step-Fourier. 

Le deuxième chapitre a été consacré à une étude générale sur les réseaux de Bragg fibrés, 
composants couramment utilisés pour leurs propriétés de réflectivité sélective en fréquence, 
comme senseurs ou comme filtre spectral de la lumière. Aussi la méthode analytique des 
modes couplés qui modélise la propagation dans ces structures périodiques a été présentée 
dans un contexte où les effets non-linéaires affectant la propagation ont été négligés. 

Le troisième chapitre présentera en premier un rappel historique sur le concept soliton. La 
dynamique de propagation des différents types de solitons optiques pouvant se propager dans 
une fibre optique et modélisée par l’équation de Schrôdinger non-linéaire comportant les 
effets linéaires et non-linéaires sera largement développée. Une autre partie du chapitre 
concernera la propagation des solitons de Bragg dans la structure FBG. Aussi, l’utilisation de 
la méthode numérique Split-Step-Fourier Method (SSFM) dans notre contribution concernera 
l’étude de l’impact des différents effets qui affectent la transmission optique par soliton. 

Le dernier chapitre a été totalement consacré à la présentation de nos différents résultats 
obtenus dans ce travail et qui concerne l’étude de la dynamique non-linéaire de la 
propagation des ondes solitaires dans un réseau de Bragg fibré dans un milieu à non linéarité 
élevée. 

Finalement une conclusion détaillée et des perspectives viendront clôturer ce travail. 
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1.1. INTRODUCTION 

Depuis les années 1970 et avec l’invention du laser en 1960, la fibre optique a attiré de plus 
en plus d'attention, en raison de son importance scientifique et technologique. Dans le 
domaine de la communication, elle représente assurément le meilleur moyen actuel pour 
transporter de très hauts débits d'informations numériques. Cette technologie constitue 
aujourd’hui une alternative rentable pour les transmissions de données sur de longues 
distances, grâce à de nombreux avantages (bande passante très grande, pertes très faibles, 
immunité accrue en regard des interférences électromagnétiques, avantages en matière de 
confidentialité et de sécurité, etc.) [1,2]. En plus de son rôle de guidage, La fibre optique est à 
la base de la réalisation de plusieurs composants optiques, tels que les filtres, les réseaux de 
Bragg, les amplificateurs, les capteurs, les modulateurs, les coupleurs,... etc. 

En effet, la propagation d’une impulsion lumineuse de forte intensité dans une fibre optique, 
peut générer un grand nombre d’effets non linéaires dont la diffusion Brillouin stimulée 
(SBS), la diffusion Raman stimulée (SRS) et les effets Kerr optiques. L’effet Kerr optique, 
qui est essentiellement une modification de l’indice de réfraction dépendant de la puissance 
optique, mène à de nombreux effets secondaires, tels que l’auto modulation de phase (SPM) 
et l’intermodulation de phase (XPM), le mélange à quatre ondes (FWM) ou l’instabilité de 
modulation (MI) [3].Ces différents effets non-linéaires compliquent énormément l’utilisation 
de la fibre mais d’un autre côté ils sont à l’origine de plusieurs applications intéressantes, 
telles que la régénération tout-optique, la conversion longueur d’onde, la génération de super 
continuum [3 ,4], etc. Cependant l’utilisation des dispositifs à base de fibres optiques non- 
linéaires a révolutionné les systèmes de communications et les besoins dans ce domaine vont 
probablement encore augmenter très fortement dans un avenir proche. 

Dans ce chapitre nous donnons un rappel général sur les propriétés fondamentales, physiques 
et géométriques des fibres optiques. En premier, nous allons présenter les différents 
paramètres caractérisant la transmission par fibre optique telles que l’atténuation et la 
dispersion puis la deuxième partie de ce chapitre sera consacré à l’étude théorique de la 
propagation d’une impulsion lumineuse dans une fibre optique décrite par l’équation de 
Schrôdinger non linéaire (NLSE). Apres présentation des notions de base sur l’optique non 





Chapitre : I 


Généralités sur l’optique non linéaire 


linéaire nous détaillerons particulièrement les principaux effets non linéaires ayant lieu dans 
la fibre optique. 


1.2. STRUCTURE DE BASE D’UNE FIBRE OPTIQUE 

La fibre optique est un guide d’onde diélectrique de forme cylindrique, composé d’un cœur de 
silice (SiOi), entouré d’une gaine également en silice. L’ensemble est généralement recouvert 
par une gaine protectrice en plastique, figure (1.1). 



Figure (1.1). Structure d’une fibre optique [5]. 

Le guidage de la lumière dans une fibre optique est assuré par une succession de réflexions 
totales des rayons lumineux à l’interface cœur-gaine, condition qui sera respectée si l’indice 
de réfraction du cœur est légèrement plus élevé que celui de la gaine (différence de quelques 
millièmes). Le dopage du cœur avec l’oxyde de germanium Ge02 ou l’oxyde de phosphore 
P 2 O 5 permet l’augmentation de son indice de réfraction alors que l’ajout de la Fluorine pour 
la gaine permet de le diminuer. Il y a plusieurs types de fibres optiques qui différent les unes 
des autres selon leur géométrie et leurs propriétés optiques et les plus utilisées sont les fibres 
monomodes et les fibres multimodes à saut ou à gradient d’indice. Pour les fibres « à saut 
d’indice », l’indice est constant et homogène dans le cœur et la gaine mais varie de manière 
discontinue à leur interface. Dans le cas des fibres « à gradient d’indice », l’indice décroit de 
manière continue du centre jusqu’à l’interface cœur-gaine [5]. 

La figure (1.2) montre la propagation de la lumière dans les trois différents types de fibres 
avec l’évolution de leurs profils d’indices de réfraction. On constate que la fibre monomode 
possède un cœur très fin (~10pm) et n’admet ainsi qu’un seul mode de propagation et permet 


y 
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une bonne transmission dans la fibre. En raison de leur faible atténuation elle reste la plus 
utilisée à l'heure actuelle pour des transmissions de données rapides et sur de grandes 
distances. Par contre les fibres multimodes (à saut ou à gradient d’indice) ont un cœur de 
diamètre plus important (50 pm -200 pm), qui permet la propagation de plusieurs modes avec 
une atténuation importante et leur utilisation reste limitée pour des bas débits ou des courtes 
distances. 



Figure (1.2). Différents types de fibre optique (a) fibre multi modes à saut d’indice, 
(b) fibre multi modes à gradient d’indice, (c) fibre monomode [5]. 


1.3. PRINCIPALES CARACTERISTIQUES D’UNE FIBRE OPTIQUE 
1.3.1. Atténuation 

Au cours de sa propagation dans la fibre optique, la puissance lumineuse d’une impulsion 
s’atténue progressivement. Cette atténuation s’exprime par une valeur en (dB/km).Soient P E 
et les puissances à l’entrée et à la sortie d’une fibre de longueur L. L’atténuation linéaire se 
traduit par une décroissance exponentielle de la puissance en fonction de la longueur de fibre : 


P s =P E e aL 


( 1 . 1 ) 
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Où a représente le coefficient d’atténuation linéique (m ’). On utilise souvent le coefficient 
adB exprimé en (dB/km) et relié à a par : 


10 i 

a d B = -— lo ë 


é p A 

1 s 

P 

V ‘E 


= 4.343a 


( 1 . 2 ) 


L’atténuation de la fibre a plusieurs origines [3-6]. On citera notamment l’absorption 
intrinsèque de la silice, la diffusion Rayleigh, la diffusion par les inhomogénéités de la fibre et 
l’absorption par des impuretés présentes dans la fibre, majoritairement les ions OH”. 
L’atténuation a également une origine mécanique : courbures, micro-courbures ou 
fluctuations du diamètre de la fibre mais aussi les épissures (zones de jonctions entre la fibre 
et un autre élément de la liaison donnant lieu à de la diffusion ou de la réflexion). 


La figure (1.3) présente la courbe d’atténuation d’une fibre en silice standard monomode. 



Figure (1.3). Courbe d’atténuation dans une fibre standard en silice monomode [6]. 

Comme on peut le voir sur la figure (1.3), l’atténuation varie suivant la longueur d’onde. La 
diffusion Rayleigh limite les performances dans le domaine des courtes longueurs d’onde 
(domaine du visible et du proche infrarouge) et le pic d’absorption, dû à la présence d’ions 
OH dans la silice, limite les performances autour de 1400 nm. Dans ces conditions 
l’atténuation typique des fibres classiquement utilisées en télécommunications est de 0.2 
dB/km à 1550 nm [3, 5, 7]. 


y 


10 








Chapitre : I 


Généralités sur l’optique non linéaire 


1.3.2. Dispersion chromatique 

Au cours de sa propagation dans une fibre optique, l’impulsion lumineuse subit un 
phénomène de dispersion qui se traduit par un étalement temporel de celle-ci, figure (1.3). 



Figure (1.3). Effet de dispersion sur l’impulsion optique. 


Plusieurs types de dispersion existent, contribuant tous à l’étalement de l’impulsion au cours 
de sa propagation dans le guide : on citera la dispersion de polarisation, la dispersion 
intermodale et la dispersion chromatique d’où cette dernière est la somme de deux effets (la 
dispersion du matériau qui dépend du type de silice, et la dispersion du guide qui dépend de la 
géométrie du guide). Par rapport aux autres formes de dispersion, c’est la dispersion 
chromatique qui joue un rôle prédominant dans les fibres optiques monomodes dont le profil 
d’indice de réfraction est de révolution asymétrique. 

Théoriquement la dispersion chromatique se traduit par la dépendance de l’indice de 
réfraction n(co) à la longueur d’onde. Chaque composante spectrale se propage à sa propre 
vitesse v=c/ n(co), ce qui induit un étalement temporel de l’impulsion au cours de sa 
propagation. La dispersion chromatique est exprimé en Ps/ (nm.km), et son effet peut être 
évalué par le développement en série de Taylor de la constante de propagation autour de la 
pulsation centrale de l’impulsion co 0 comme suit [3] : 

/3{co) = n(co).col c (1.3) 


: A) + /?, (co — CO Q ) + 


R (r 


Avec : 


y 


ii 





Chapitre : I 


Généralités sur l’optique non linéaire 


Pm 


d m p 

dco’" 


CO=G)q 


m = 0,1, 2,3,.... 


(1.5) 


Les paramètres /?, ,J3 2 représentent respectivement l’inverse de la vitesse de groupe v g et la 
dispersion de la vitesse de groupe (GVD). Les deux sont reliés à l’indice optique et ses 
dérivées par : 


A 


A 




i dv g 

v 2 g dco 


C 


dn 

dco 


co 


dn 
dco 2 




( 1 . 6 ) 

(1.7) 


Généralement en télécommunication optique on parle du coefficient de dispersion D qui est 
relié avec le paramètre de dispersion de deuxième ordre par : 

dU)=- 2^Æ d-8) 

Pour une longueur d’onde particulière A 0 . la courbe de dispersion s’annule. Cette longueur 
d’onde de dispersion nulle permet de définir deux régimes de dispersion [3]. Le premier 
régime dit nonnal c’est pour les longueurs d’ondes X <A 0 (|L>0). Dans ce régime les 
composantes spectrales de haute fréquence se propagent plus lentement que celles de basses 
fréquences. C’est le contraire qui se passe dans le deuxième régime dit anormal pour les 
longueurs d’ondes X >X <) (p2<0). 


1.4. OPTIQUE NON LINEAIRE 
1.4.1. Introduction 

La compréhension des effets non linéaires dans les fibres optiques est basée sur l’étude 
théorique de la propagation d’onde optique dans un milieu non linéaire et dispersif. Dans la 
théorie optique, fonde lumineuse est une onde électromagnétique composée d’un champ 
électrique et d’un champ magnétique, dans le cadre de cette théorie, l’évolution du champ 
optique dans un milieu matériel est gouvernée par les équations de Maxwell [3,8] : 

V ■ D = p ( L9 ) 


12 




Chapitre : I 


Généralités sur l’optique non linéaire 


V x H = J + — 
dt 

Y7 P ^ 

Vx E =- 

dt 

VB = 0 


( 1 . 10 ) 


( 1 . 11 ) 


( 1 . 12 ) 


J et p sont la densité de courant et la densité de charge du milieu. E et H sont respectivement 
les champs électrique et magnétique, alors que D et B sont l’induction électrique et 
l’induction magnétique. On a aussi les équations suivantes : 


D = s 0 E + P 
B = jU 0 (m + h) 


(1.13) 


(1.14) 


Où P et M sont respectivement la polarisation et la magnétisation. La permittivité diélectrique 
du vide (£o =8.854.10 12 F/m), et la perméabilité magnétique du vide (po= 4îi. 10 7 H/m) sont 
reliés à la vitesse de lumière dans le vide (c=299 792 458 m/s) par cette relation : 


Mo e 0 — 


(1.15) 


Pour les matériaux non magnétiques (M =0) comme dans notre cas (fibres optiques en silice) 
la densité de courant et la densité de charge sont nulle (J = 0, p = 0) et avec la combinaison 
des équations (1.9-14) et l’utilisation de la relation (1.15), on obtient : 


VxVx£(?,/)=-^- 


1 d 2 E(r,t) d 2 P(r,t) 

2 ^ 2 MO 2 

: dt dt 


(1.16) 


VxVx£=V 


.(vJ) 


-v 2 e 


(1.17) 


On substituant la relation (1.17) dans l’équation (1.16), on obtient la forme de l’équation 
d’onde qui décrit la propagation de la lumière dans les fibres optiques [9]. 


u2îr „ , 1 d 2 E(r,t) d 2 P(r,t) 

~ VE(r - = 


(1.18) 
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En effet, le procédé le plus commun pour décrire les phénomènes optiques non-linéaires est 
basé sur l’expression de la polarisation induite p en termes de force du champ électrique 
appliqué E . La raison pour laquelle la polarisation joue un rôle clef est qu’elle peut agir 
comme source de nouvelles composantes du champ électromagnétique. 


1.4.2. Polarisation linéaire et non linéaire d’un milieu 

Un milieu matériel peut être considéré comme un ensemble de particules chargées, des 
électrons de charge négative et des ions de charge positive. Si le mi lieu est conducteur, 
l’application d’un champ électrique entrainera un mouvement des charges donnant lieu à un 
courant électrique. Dans le cas d’un matériau diélectrique, la liaison entre les charges étant 
très forte, seul un faible déplacement est permis et ces petits mouvements de charges 
(électrons dans un sens et ions dans l’autre sens) induisent dans le matériau des moments 
dipolaires électriques .c.-à-d. l’effet d’un champ électrique sur un diélectrique est d’induire 
une polarisation, sous la forme P — f{E ). Pour de faibles intensités du champ appliqué 
devant le champ atomique du milieu et suffisamment éloigné de ses fréquences de résonances, 
ces dipôles induits sont proportionnels au champ électromagnétique. Nous sommes alors en 
présence d’une polarisation linéaire. Cependant, pour un champ intense (qui dépasse 
typiquement le seuil 10 6 V/m, ceci n’est vérifié que pour des sources laser), les dipôles ne 
croissent plus linéairement avec le champ appliqué, mais réagissent de manière non-linéaire. 
La polarisation induite par les dipôles électriques est une fonction non-linéaire du champ 
électrique, qui fait intervenir la notion de susceptibilité optique non-linéaire, donnée par 
l’équation suivante [3,9] : 

P = s 0 (Z m -E + x ay ÉÉ + x 0) ':ÉÉÉ + .) (L19) 

Où j tn) est la susceptibilité d’ordre n, (n=l, 2,3,..). La polarisation comporte deux termes : 

P = P L + P NL 
? L =e 0 Z m -E 

p NL =8 0 z i2y ÉÊ + s 0 z O) \ÉÉÊ + 


( 1 . 20 ) 

( 1 . 21 ) 

( 1 . 22 ) 
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/S : est la contribution linéaire à la polarisation, /^exprime la réponse non linéaire du milieu. 
Pour une fibre optique en silice (S KL) tous les termes pairs de la susceptibilité non linéaire 
sont nuis pour des raisons d’une symétrie d’inversion au niveau moléculaire. Cependant la 
polarisation dans les fibres optique prend la forme générale suivante : 

P = s 0 (Z a) ■É + x 0) :ÉÉÉ + .) ( L23 ) 

Les parties réelles et imaginaires de la susceptibilité linéaire d’ordre 1 sont liées à l’indice de 
réfraction et au coefficient d’atténuation par les relations suivantes [3] : 


n(a>) = l + ^-Re[j (1, (ru)] 

(1.24) 

a(œ) = — lm[x a \œ)] 

(1.25) 


ne 


La susceptibilité non-linéaire d’ordre 3 est à l’origine de plusieurs effets non-linéaires qu’on 
peut diviser en deux catégories [3,10]. Les effets élastiques et les effets inélastiques. Pour les 
effets élastiques il n’y a aucun échange d’énergie entre le faisceau laser et le matériau dans 
lequel il se propage (tel que l’effet Kerr). Par contre dans les effets inélastiques (nommés en 
général diffusion) il y a un échange d’énergie entre le champ optique et le matériau. 


1.4.3. Effets non- linéaires dans la fibre optique 
I.4.3.I. Effet Kerr optique 

La polarisation non linéaire induite dans la fibre optique par la propagation d’un champ 
optique de forte intensité, est responsable d’une dépendance en intensité de la valeur de 
l’indice de réfraction. Cette dépendance de l’indice de réfraction vis-à-vis de l’intensité est 
connue sous le nom d’effet Kerr optique [3,7 ,11]. On peut écrire la relation (1.23) qui 
exprime le développement du vecteur de polarisation en série de puissance du champ 
électrique sous la forme suivante: 

P = s 0 { X m +s NL )Ë ( L26 ) 



Avec 
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£ nl 



(1.27) 


La polarisation et l’indice de réfraction sont étroitement liés par la relation de Maxwell (1.13) 
comme suit : 


D = s Q Ê + P = n 2 Ë 


(1.28) 


Cette relation conduit à : 

n 2 =l + z n) + s NL = K + A n NL ) 2 « n 2 + 2n 0 An NL (1.29) 

Donc, l’indice de réfraction est composé d’une partie linéaire no, et d’une partie non linéaire 
Àn N L, qui sont données par : 


i/ 
(1> V2 


no=(X + Z w ) 

s 


A n NL ~ 


'NL 

2 n n 


Il vient finalement, qu’à une pulsation co donnée on a [3,12] : 


n(a>, E ~) = n Q (co) + n 2 E 


(1.30) 

(1.31) 


(1.32) 


Où ri 2 représente le coefficient d’indice de réfraction non linéaire donné par : 


n 2 = 


8 n, 


X 


( 3 ) 


(1.33) 


L'indice de réfraction peut également s'exprimer en fonction de l'intensité optique I par : 


n(co,I) = n 0 (û)) + n 2 I (1-34) 

Où I — a\E\ , avec a = e 0 cn 0 / 2, et n!, = /a , n^ est appelé l’indice Kerr. L’effet Kerr 

optique, ou la dépendance de l’indice de réfraction vis-à-vis de l’intensité est l’effet le plus 
important dans la fibre optique. Cet effet est à l’origine d’autre effets non linéaires 
intéressants que nous détaillerons par la suite [1-7], comme l’auto-modulation de phase, la 
modulation de phase croisée, le mélange à quatre ondes, l’instabilité modulationnelle et 
encore l’existence d’impulsions solitoniques. 
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I.4.3.1.a. Auto-modulation de phase (SPM) 

La modification non linéaire de l’indice de réfraction sous l’effet d’un champ optique de forte 
intensité, entraine un déphasage non linéaire auto-induit sur l’onde se propageant dans la 
fibre. Cet effet est nommé auto-modulation de phase (SPM), donnée par [3 ,13]: 

^Ûspm = n 2 k 0 IL (1.35) 

Avec k 0 le vecteur d’onde, I l’intensité lumineuse, L la longueur de la fibre. L’effet d’auto¬ 
modulation de phase conduit en général à un élargissement spectrale se traduisant par la 
dépendance temporelle du décalage fréquentiel appelé (chirp) défini par : 

d^spM (T) (1.36) 

dT 



I.4.3.1.b. Modulation de phase croisée (XPM) 

On appelle la modulation de phase croisée le déphasage non linéaire supplémentaire induit 
par l’intensité d’autres ondes de longueur d’onde ou de polarisation différentes, se propageant 
simultanément avec la première onde. Ce déphasage est donné par [3,14] : 

A</> xpm = 2n 2 k 0 I 2 L (1.37) 

avec h, l’intensité de la deuxième onde. 


I.4.3.I.C. Mélange à quatre ondes (FWM) 

La modulation d’indice de réfraction via l’effet Kerr optique induit un phénomène physique 
de transfert d’énergie entre les fréquences qui se propagent dans la fibre optique et qui 
donnent heu à l’apparition de nouvelles fréquences. Cet effet dit le mélange à quatre ondes 
(FWM), peut être décrit par l’annihilation de deux photons (dits de pompe) de même 
fréquence (processus dégénéré) ou de fréquences déférentes wq et w 2 (processus non 
dégénéré), et la création simultanée de deux autres photons de fréquences caractéristiques w 5 
et W4 tel que wq + w 2 = W3 + w 4 . L’efficacité de FWM est liée au respect d’une condition 
dite d’accord de phase entre les différents vecteurs d’ondes mis enjeu. Cet effet est à la base 
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d’autres effets non linéaires comme l’amplification paramétrique et l’instabilité 
modulationnelle [15]. 


I.4.3.1.d. Instabilité modulationnelle (IM) 

L’instabilité modulationnelle (IM) est un phénomène non linéaire qui apparait lorsqu’une 
onde se propage dans un milieu non linéaire et en régime de dispersion anormale c.-à-d. le 
coefficient de dispersion B 2 est négatif. Cet effet se traduit par la modulation en amplitude de 
l’onde sous l’effet d’une perturbation. Lorsque l’onde incidente est initialement modulée, 
l’instabilité peut amplifier cette modulation, on parle alors d’instabilité induite. Si l’onde 
incidente est continue et si la modulation croit à partir de bruit, on parle d’instabilité 
spontanée [16]. 


I.4.3.2. Diffusions inélastiques Raman et Brillouin stimulées 

Les diffusions Raman (SRS) et Brillouin (SBS) stimulées, sont des effets non linéaires qui 
apparaissent lorsque la fibre subit un champ optique de forte intensité. Ces deux phénomènes 
consistent à l’excitation des modes de vibration de la silice aux niveaux moléculaires. Dans le 
cas de la diffusion Raman, une interaction se fait entre l’onde lumineuse et les phonons 
optiques de la silice. Mais la diffusion Brillouin résulte de l’interaction entre l’onde optique 
et les phonons acoustiques, donnant lieu à l’apparition d’une onde Stokes décalée en 
fréquence par rapport à l’onde incidente (il est environs de 13 THz pour l’effet Raman, et 
environs de 11GHz pour l’effet Brillouin) [17]. La valeur de ce décalage varie d’une fibre à 
une autre selon le dopage utilisé. On parle d’onde anti-Stokes pour la génération d’une onde 
décalé vers les hautes fréquences comme le montre la figure (1.4). 

On note que l’effet Raman est un effet très important, il est à l’origine des amplificateurs à 
large bande spectrale et des lasers accordables en longueur d’onde. 
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Figure (1.4). La génération d’onde Stokes et anti-Stokes par l’effet Raman [17]. 


1.5. EQUATION DE PROPAGATION NON-LINEAIRE 

La propagation d’une impulsion lumineuse intense dans une fibre optique se modélise par 
l’équation de Schrôdinger non linéaire (NLSE). Elle est obtenue à partir des équations de 
Maxwell et de l’équation d’onde (1.18), avec l'inclusion des termes de polarisation non 
linéaire qui sont responsable de plusieurs effets non linéaires [3]. L’équation de Schrôdinger 
peut décrire tous les effets linéaires et non linéaires qui affectent l'impulsion au cours de sa 
propagation dans la fibre optique. 

Pour dériver l’équation de Schrôdinger on suppose la fibre monomode, idéalement 
homogène et isotrope, caractérisée par une faible non-linéarité et une forte dispersion (i.e. Pnl 
«< Pl). Aussi on suppose que la largeur spectrale de l’impulsion étudié est très inférieure à 
sa fréquence centrale (Àco/coo << 1 ), pour des impulsions courtes (7o<0.1 ps). Dans ces 
conditions l’approximation de l’enveloppe lentement variable est valide et la non-linéarité 
peut être traitée comme phénomène instantané [3]. 

En considérant la propagation des impulsions suivant une seule dimension z, l’équation 
(1.18) peut se réécrire dans le domaine fréquentiel comme suit : 

~ a 2 co 2 
_dz 2 + c 2 

Par l’utilisation de la transformée de Fourier suivante : 


y 
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avec : 


+00 

Ë(z,œ)= | E(z,t)e i<a dt 


—00 

E(z,t) = A(z,t)e- ia » 


(1.39) 

(1.40) 


Dont l’enveloppe A(z, t) varie lentement dans l’espace et dans le temps devant la 
porteuse œ 0 . En effet on peut séparer le champ électrique E(x, y, z, ùj) en deux parties, partie 
longitudinale et une autre transversale tel que : 


Ë(x,y l z,<D ) = F(x, y, n>)Ë(z, n>) 


(1.41) 


où la partie transversale est solution de l’équation de Helmholtz suivante [3,18]: 


ô 2 | a 2 

ôx 2 + ôy 2 


+ k~(oo) 


P{œ) 2 F 


(1.42) 


/?(n>) : est la constante de propagation développée en série de Taylor dans l’équation 
(1.4). £(u>) : est la fonction diélectrique donnée par cette relation : £(n>) — 1 + X 1 (m). 

La fonction diélectrique liée à l’indice de réfraction et au coefficient d’atténuation vient 
comme suit : 


s(co ) =(n + iac/ 2 cof 

k{co) = tJs(co)û)/ c = n L co/c + ia/2 = f3(co) + ia/2 


(1.43) 

(1.44) 


En introduisant k(oS) , l’équation (1.38) prend la forme suivante : 

d 2 


dz‘ 


+ k (aj) 


E = —/JqCD Pj 


NL 


(1.45) 


L’utilisation de l’équation (1.40) pennet de développer l’équation (1.45) en une équation de 
premier ordre comme suit : 


d 


- ik(ùj) 

dz 


A = iju {) or P NL / 2 k{œ)A eff 
Où A e ff est l’aire effective du mode guidé [3]. 


(1.46) 


11| F(x,y,co) dxdy 


A — v 

'X# +00 


(1.47) 


1 F<* , y,<y)| dxdy 





Chapitre : I 


Généralités sur l’optique non linéaire 


Les polarisations linéaire et non linéaire du milieu sont données par ces relations [3] : 

P L (z,t) = s 0 \ R [X) (t)E(zJ - T)dr (!- 48 ) 

J—00 

P NL (z,t)=£ 0 E(z,t) f R {3) (r)\E(z,t- t)\~ dr (1.49) 

J—CO 1 1 

Où r — t — t' , R 1 et R 3 sont les fonctions réponse linéaire et réponse non 
linéaire successivement. La transformée de Fourier de R 1 (t) est la susceptibilité 
linéaire X 1 (m), et la transformée de Fourier de R 3 (t) est la susceptibilité non linéaire 
d’ordre trois X 3 (u>), dont la fonction de réponse non linéaire R 3 est définie par cette 
relation [3,19]: 

R°< (r) = 2n 2 n L [a-f R )<5(r) + f R h R (r)] < 150 ) 

Cette fonction contient la non-linéarité de Kerr instantanée 8 (r), la fonction réponse Raman 
retardée h R ( T ), et la contribution de l’effet Raman à l’effet Kerr optique f R . En introduisant 
les relations (1.44), (1.49) et (1.50) dans l’équation (1.46) et avec quelques simplifications 

mathématiques et l’utilisation de la transformation suivante :T = (t — = (t — B x z) , on 

obtient dans le domaine temporel l’équation de propagation non linéaire suivante [3] : 


d . ~ a 
- iD H- 

ôz 2 


| A(z, T ) = iy(oj)(\ + — A) x A(z, T) 
a> Q oT 




(1.51) 


Où y (en) est le coefficient non linéaire défini par cette relation [3] : 


Y(œ) = 


n 2 co 


cA 


(1.52) 


eff 


L’opérateur D est appelé opérateur dispersion, donné par cette relation 


^ ,-m+l a mo am 

D = (.j~ïi)o>=o) 0 


(1.53) 


L’équation (1.51) est nommée l’équation non linéaire de Schrôdinger généralisée. Elle décrit 
dans le domaine temporel la propagation de l’enveloppe A(z,t) de l’impulsion optique 
suivant l’axe longitudinal z de la fibre optique. Cette équation prend en compte les effets 
linéaires présentés par le paramètre de perte (a), l’operateur (D) de dispersion, les effets non 
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linéaires de type Kerr décrits par le coefficient (y) , les diffusions inélastiques décrites par la 
fonction h R(t j et le terme proportionnel à coo~ l responsable de l’effet d’auto-raidissement. 


1.5.1. Equation de Schrodinger non linéaire simple 


Si on néglige les effets Raman et Brillouin et les ordres supérieurs de dispersion chromatique 
(pour des impulsions longues, T 0 > 5 ps), l’équation (1.51) de Schrôdinger non linéaire 
généralisée prend la forme simple suivante [3,20] : 


dA a .B^ d 2 A . , ,2 

— + — A + i — —--iyA A = 0 
dz 2 2 dT 2 ' 


(1.54) 


Cette équation est à la base de nombreux travaux de recherche dans le domaine des 
télécommunications, elle a été démontrée dans les fibres optiques en 1973 par Hasegawa 
et Tappert [3]. Elle permet de décrire des phénomènes aussi variés que la dispersion 
chromatique d’ordre 2, l’instabilité modulationelle, l’auto-modulation de phase ou encore la 
propagation de soliton. 


On peut écrire l’ENLS simple en d’autre forme où les paramètres sont normalisés et les pertes 
négligées [3]. Cette forme normalisée est très utile dans le traitement de la transmission par 
soliton: 


. du 
1 — 
dZ 


B 2 1 d"u 1 12 

■ + \u\ u = 0 


B 2 2 dx t 


(1.55) 


Où t — T/T 0 est le temps normalisé, To est la largeur v à mi-hauteur de l’impulsion, Z est la 
distance normalisée donnée par : Z = z/L D , L D est la longueur de dispersion à partir de 
laquelle les effets de dispersion commencent à devenir importants. Elle est définie par cette 
relation [3]: 


T 

D “ I B, 


(1.56) 


Avec u l’amplitude normalisée définie par : 

u = N — 


V / l d a 


(1.57) 
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Où Po est la puissance crête de l’impulsion A (z, t), N est l’ordre du soliton défini par : N 2 — 
L d /L nl où L nl est la longueur non linéaire à partir de laquelle les effets non linéaires 
deviennent significatifs. Elle est donnée par : 


-‘NL 



(1.58) 


1.5.2. Différents régimes de propagation 

Selon les valeurs de L nl et L D et la longueur de la fibre L, on distingue quatre régimes de 
propagation [3]: 

1. Régime linéaire non dispersif : pour ( L nl » L,L d » V) , les effets dispersifs ou 
non linéaires ne jouent aucun rôle significatif, l’impulsion conserve sa forme durant la 
propagation. Dans ce cas l’équation (1.54) prend la forme suivante : 



& 


2. Régime dispersif : pour (L nl » L,L d < L), dans ce cas les effets non linéaires sont 
négligeables, l’impulsion subit un élargissement temporel sous l’effet de la dispersion. 
L’équation (1.54) devient : 

. du 1 d 2 u 

i — +- T = 0 

ôz 2 dr 2 

3. Régime non linéaire: lorsque (L D » L,L nl < L) , la dispersion est négligeable, 
l’impulsion subit un changement de phase et l’équation (1.54) devient : 

■ 2 

u\ u = 0 

4. Régime soliton : lorsque ( L D S L nl < L), dans ce cas les effets non linéaires et l’effet 
de dispersion agissent simultanément et se neutralisent. Les impulsions restent stables 
sans déformation au cours de la propagation, ceux sont des ondes solitaires. Dans ce 
cas la propagation est modélisée par l’ENLS (1.55). 


. du 
l -b 

dz 
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1.5.3. Différents types d’impulsions optiques 

L’enveloppe A(z, T) de l’impulsion optique injectée dans la fibre prend en général les 
différentes formes suivantes [3]: 

, - T 2 

a) impulsion gaussienne : le profil de l’impulsion est de type :^4(0, T) = V P 0 (—-) où 

T 0 est la demi-largeur (au point d’intensité P 0 /e), elle est reliée avec la largeur totale 
à mi-hauteur par la relation T FWHM — 2(/n2) 1 / 2 7’ 0 « 1 ,66ST 0 . Si l’impulsion initiale 
est chirpée linéairement, son profil se présente comme 

suit: A(0,T) = tJP o ( -—,— où C est le paramètre du chirp (dépendance 

temporelle de la fréquence instantanée). 

b) impulsion de type sécante hyperbolique : l’enveloppe du champ incident prend 
cette forme : ,4(0,7) = J~P^sech(?-)exp(- , avec T FWHM = 2 ln(l + VI) T 0 « 

'0 ^'o 

1.7637 0 . La forme de type sécante hyperbolique est considérée en particulier dans le 
contexte des solitons optiques et des impulsions émises par des lasers à modes 
bloqués. 


1.6. RESOLUTION NUMERIQUE DE L’EQUATION DE SCHRÔDINGER 

En effet l’équation de Schrôdinger non linéaire est une équation différentielle partielle non 
linéaire, où elle est très difficile à résoudre analytiquement. Une approche numérique est donc 
souvent nécessaire pour mieux analyser l’impact des effets non linéaires sur la propagation 
des impulsions dans la fibre optique. Parmi un grand nombre de méthodes numériques 
utilisées on trouve la méthode Split-Step-Fourier (SSF). Elle a été utilisée pour la première 
fois par Hasegawa en 1973, puis elle a été employée intensivement du fait de son efficacité et 
de sa précision dans les calculs [3,21]. Dans ce mémoire, une grande partie des simulations 
numériques présentées ont été réalisées à l’aide de cette méthode et grâce à des algorithmes 
développés sous un environnement Matlab. 
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Le principe de cette méthode est de considérer l’équation de Schrôdinger non linéaire comme 
pouvant se mettre sous la forme suivante [3 ,22]: 


PA 

i — - (D + N)A 

dz 


(1.59) 


Où D et N représentent l’opérateur linéaire et l’opérateur non linéaire. Si on prend comme 
exemple l’équation simple de Schrôdinger (1.54), les deux opérateurs s’écrivent comme suit : 


1 d 2 
D = —(a + iB 7 —-) 

2 dT 2 


N =iy\A\ 2 


(1.60) 

(1.61) 


Généralement la dispersion et la non linéarité agissent simultanément suivant la longueur de 
la fibre. Cependant la méthode SSF est basée sur une approximation qui consiste à dire que 
sur de très courtes distances de propagation h, les operateurs D et N commutent. Les effets 
linéaires et non linéaires agissent indépendamment l’un de l’autre comme le montre la figure 
(1.5). 



Figure (1.5) : Représentation schématique de la méthode Split-Step-Fourier. 


La propagation entre les points z et z+ h se fait donc en deux étapes : 

1. seuls les effets linéaires agissent (D ^ 0 et N — 0). 

2. seuls les effets non linéaires agissent ( N ^ 0 et D — 0). 

Dans ces conditions, la résolution de l’équation (1.54) conduit à deux étapes : 


y 
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1) . Etape dispersif : 

On a : (D A 0 et N — 0). L’équation (1.54) est devient : / ^ = f)A(z,t) 

ôz 

La résolution de cette équation dans le domaine fréquentiel se présente comme suit : 

Â(z + h, co) = Â(z, co) exp (hD) (1-62) 

Comme l’opérateur différentiel ^ est équivalent à une multiplication par - ico dans le 
domaine fréquentiel, donc l'expression fréquentielle de l’opérateur linéaire D est donnée par : 

D =- — {a - iB 2 co 2 ) e t avec l’utilisation de la transformé de Fourier inverse TF 1 , on peut 
écrire la solution (1.61) dans le domaine temporel comme suit : 

A'(z + h,T)= TF 1 { À{z,co) exp(-^ (a — iB 2 co 2 )) } (1.63) 

2) . Etape non linéaire : 

L’opérateur N s’applique ensuite dans le domaine temporel : 

A(z + h,T) = A , (z,T)exp(hN) C 1 - 64 ) 

Au final, la solution de l’équation de Schrôdinger (1.53) devient : 

A(z + h,T) = exp(ih/\A\~)TF^ 

En répétant ces deux opérations sur toute la longueur de la fibre, nous réalisons ainsi une 
simulation numérique de la propagation des impulsions optique dans une fibre optique. 


h 

2 1 


A(z, co) exp(- — (a - iB 2 co~ )) 


(1.65) 


1.7. CONCLUSION 

Dans ce chapitre nous avons étudié théoriquement la propagation des impulsions lumineuses 
dans la fibre optique comme guide d’onde dispersif et non linéaire. Cette propagation peut 
être modélisée par l’équation de Schrodinger non linéaire dérivée à partir des équations de 
Maxwell. La résolution de cette équation étant très difficile analytiquement d’où une 
résolution numérique par la méthode Split-Step Fourier s’avère efficace et donne des 
résultats rapides et précis. 
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1.1. INTRODUCTION 

Les fibres optiques peuvent changer leurs propriétés optiques quand elles sont exposées à un 
champ optique de forte intensité (UV) et ce phénomène a été découvert pour la première fois 
par Hill et al. en 1978 [1]. La photosensibilité de la silice induit une modulation périodique de 
l’indice de réfraction le long de la fibre et cette périodicité se traduit par la formation d’un 
réseau de Bragg fibré (FBG), au cœur de la fibre. En 1989 l’équipe de G. Meltz et W. Morey, 
a pu décrire pour la première fois une méthode de photo-inscription des réseaux de Bragg 
dans les fibres optiques [2]. En l’espace de quelques années, la technologie des réseaux de 
Bragg photo-inscrit a révolutionné le domaine des transmissions optiques et le domaine des 
capteurs, grâce à plusieurs applications. 

Ce chapitre sera totalement consacré à l’étude des réseaux de Bragg fibrés. Les principes 
physiques, les différentes méthodes d’inscriptions ainsi que les différents domaines 
d’applications seront présentés et sera particulièrement détaillée la théorie des modes couplés 
qui décrit la modélisation de la propagation des ondes lumineuses dans les réseaux de Bragg. 


II.2. DEFINITION D’UN RESEAU DE BRAGG 

La forme simple d’un réseau de Bragg fibré (FBG), consiste en une modulation périodique de 
l’indice de réfraction du cœur d’une fibre monomode selon l’axe de la fibre (Figure II. 1) [3]. 


a 

o 

"0 

ai 

c 





Figure (II.l) : Principe d’un réseau de Bragg [3]. 
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En effet le (FBG) se comporte comme un filtre de longueur d’ondes. Pour les longueurs 
d’ondes situées autour de la longueur d’onde de Bragg À B , elles sont partiellement réfléchies 
par le réseau, par contre les autres elles seront transmises. Les ondes réfléchies vérifient la 
relation suivante : 


À B — 2n e ffA (II-1 ) 

Où n e ff est l’indice de réfraction effectif du mode fondamental et a représente le pas du 
réseau désignant la période de modulation de l’indice de réfraction. 


II.3. METHODES D’INSCRIPTION 

Le réseau de Bragg est une micro structure de quelques millimètres, qui peut être inscrit sur le 
cœur d’une fibre optique photosensible. C’est en injectant dans une fibre monomode avec un 
cœur dopé à l’oxyde de germanium une onde lumineuse issue d’un laser à argon ionisé (f= 
488 nm) que l’équipe de Hill et al. a réussi par hasard la première inscription d’un (FBG) en 
1978[ 1]. Le résultat a été la création d’une onde stationnaire entre les deux extrémités de la 
fibre et ainsi une variation périodique de l’indice de réfraction apparait sous l’effet de 
l’interférence entre fonde injectée et fonde réfléchie. Cette méthode a permis l’inscription 
d’un (FBG) longitudinal, très long avec une variation périodique de l’indice très faible. 
Depuis, plusieurs études de procédés d'inscription de réseaux de Bragg ont suivi cette 
découverte et la première méthode d’inscription efficace a été celle du procédé holographique 
transverse proposé par l’équipe de G. Meltz et W. Morey en 1989 [2]. Quelques années plus 
tard (1993), l’équipe de Hill et al. a proposé une nouvelle méthode plus évoluée en utilisant 
un élément diffractif appelé masque de phase [4]. Ces deux techniques sont brièvement 
décrites ci-dessous. 


II.3.1. Méthode holographique transverse 

Le principe de cette méthode d’inscription transverse est basé sur la division d’un champ 
ultraviolet selon deux faisceaux qui interfèrent, pour former des franges d’interférences sur la 
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partie de la fibre insolé et induisent ainsi une modulation de l’indice de réfraction. Le 
montage expérimental utilisé est décrit sur la figure (II.2). En utilisant des miroirs orientés 
vers la zone à graver et un champ d’interférences produit par deux faisceaux (UV) cohérents 
à 244 nm, les auteurs ont permis de graver un (FBG) transverse à la fibre, de longueur 4.4 
mm, avec une modulation d’indice de réfraction d’environ 3.10 5 [3,2]. 


UV plaque de compensation 



Figure (II.2) : Principe du montage de la méthode holographique transverse [2]. 

La longueur d’onde de Bragg du réseau photo inscrit par cette méthode étant lié à l’angle que 
font les deux faisceaux qui interfèrent, permet la fabrication de réseaux de Bragg à des pas 
différents et avec des caractéristiques appropriées pour des applications dans le domaine des 
télécommunications et le domaine des capteurs. La longueur d’onde de Bragg est donnée par 
cette relation : 


-, n e ff^ UV 

A B ~ — 9 

n uv sin- 


(II. 2) 


Où n e f f est l’indice de réfraction effectif du cœur de la fibre, Â uv est la longueur d’onde du 
laser, n uv est l’indice du milieu (l’air) d’où passe le faisceau laser. 
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II.3.2. Méthode du masque de phase 

La technique d’inscription des (FBG) par l’utilisation du masque de phase est la technique la 
plus courante du fait de son excellente reproductibilité et sa forte efficacité. Un masque de 
phase est une lame de silice sur laquelle a été gravé un réseau de diffraction. Ce dernier 
diffracte la lumière du laser incident dans les ordre +1 et -1, [3,4], figure (II.3). 

UV 


masque de phase 



Figure (II.3) : Principe du montage de la méthode du masque de phase [4]. 

l’interference des deux ordres de diffraction créent une modulation de l’indice du cœur avec 
un pas fixé par le pas du masque de phase tel que : 


a b — 


_ -d masque 


(II.3) 


Dans ce cas il est nécessaire d’utiliser un masque de phase différent pour chaque réseau et il 
est également possible d’inscrire un réseau de Bragg à pas variable (communément appelé 
réseau de Bragg chirpé), par l’utilisation d’un masque de phase à pas variable. 


II.4. DIFFERENTS TYPES DES RESEAUX DE BRAGG 

Plusieurs types de réseaux de Bragg peuvent être gravés sur la fibre optique et ils diffèrent les 
uns des autres selon le profil de l’indice de réfraction qui est décrit par une fonction quasi- 
sinusoïdale comme suit [3,5] : 
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An e// (z) = A n eff (z) {l + v(z)cos(y z + 0(z))} (II.4) 

Où A n e ff(z) est la modification moyenne de l’indice effectif, v(z) l’apodisation de la 
modulation d’indice, a la période de la modulation et 0(z) sa phase. 

La figure (II.4) présente les différents profils d’indice de réfraction pour les principaux types 
de réseaux de Bragg [3,5, 6]. 

a) Réseau uniforme : Le réseau uniforme est caractérisé par une période, une phase et 
une amplitude constante le long du réseau, 

b) Réseau à amplitude variable: An e ^(z) varie en fonction de z de façon gaussienne, 

c) Réseau apodisé : v(z) varie d’une façon sinusoïdale et A n e ^(z) est constant, 

d) Réseau à pas variable (chirpés) : la phase 0(z) est non nulle et le pas du réseau 
varie avec z, 

e) Réseau à saut de phase : c’est un réseau uniforme ayant une phase variant par 
endroit de façon discontinue par sauts, 

f) Réseau superstructuré : v(z)varie d’une façon périodique. 



(C) 



A n 


eff 



(O 


z 


Figure (II.4) : Principaux types de réseaux de Bragg [3]. 
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ILS. DOMAINES D’APPLICATIONS DES RESEAUX DE BRAGG 

Depuis le début des années 90, la technologie d’inscription des réseaux de Bragg dans les 
fibres optiques a connus un développement rapide, et a atteint une certaine maturité. De nos 
jours le réseau de Bragg est devenu le composant le plus important et stratégique de plusieurs 
applications dans la majorité des domaines d'applications de la fibre optique, en particulier 
dans le domaine des télécommunications et des capteurs [3,7]. 


II.5.1. Domaine des télécommunications 

Les réseaux de Bragg trouvent un intérêt majeur dans les systèmes de transmission optique à 
raison de leurs atouts suivants [3,7, 8]: 

•S compatibilité avec les fibres optiques utilisées dans les réseaux de télécommunication 
•S faibles pertes d’insertion et une très faible sensibilité à la polarisation. 

•S fabrication à faible coût de dispositifs de très grande qualité sélectifs en longueur 
d’onde. 

■f faciles à intégrer dans les systèmes de télécommunication. 

Z disposent d’une sensibilité élevée d’une très grande fiabilité sur la durée. 

Z une capacité à fournir des mesures absolues sans avoir besoin d’un 
référencement. 

Ces atouts en font également des composants très attractifs pour plusieurs applications 
notamment comme : filtres réglables à bande étroite et à large bande, filtres de longueurs 
d’ondes sélectifs multiplexeurs/démultiplexeurs dans les systèmes à multiplexage en 
longueurs d’ondes (WDM). La figure (II.5.a) présente l’exemple d’un filtre passe bande 
fabriqué par un réseau de Bragg associé avec un circulateur. La figure (II.5.b) présente un 
multiplexeur (OADM) fabriqué avec l’association d’un réseau de Bragg et de deux 
circulateurs. Les circulateurs sont des dispositifs coûteux, mais on peut les remplacer par 
d’autres dispositifs peux coûteux comme les interféromètres Mach-Zhender, qui sont utilisés 
pour réaliser des multiplexeurs (OADM), ou bien les interféromètres Fabry-pérot (FFP) pour 
fabriquer des filtres. 
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Bragg grating al 




Figure (II.5) : (a) Filtre passe bande à réseau de Bragg associé avec un circulateur ; 

(b) Multiplexeur à réseau de Bragg associé avec deux circulateurs. 

Les (FBG) peuvent être aussi utilisés comme compensateurs de la dispersion chromatique 
dans les réseaux de télécommunications à grande distance. La figure (II.6) montre le principe 
de la compensation de la dispersion chromatique en utilisant un réseau de Bragg chirpé en 
association avec un circulateur. 



Figure (II.6) : Compensateurs de la dispersion chromatique à réseau de Bragg chirpé. 

On trouve aussi les réseaux de Bragg (FBG) dans d’autres applications comme : analyseur de 
spectre, lasers spéciaux à bande étroite, correction du gain et amélioration du rendement de la 
pompe pour les amplificateurs optiques à fibre dopée à l’erbium, convertisseurs de mode à 
fibres optiques et codeurs et décodeurs dans les systèmes OCDMA. 


II.5.2. Domaine des capteurs. 

Un réseau de Bragg possède des caractéristiques uniques pour être utilisé comme capteur. En 
effet, lorsqu’un réseau de Bragg est soumis à une perturbation physique (température, 
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déformation, pression hydrostatique, champ électromagnétique...), il se dilate ou se contracte, 
ce qui modifie son pas, son indice de réfraction et par conséquent la longueur d’onde de 
réflexion ou de transmission. Cependant le principe de fonctionnement d’un capteur à réseau 
de Bragg est de mesurer le décalage de la longueur d’onde de Bragg par rapport a une 
référence choisie [5, 9,10], figure (II.7). 


Signal réfléchi 



Figure (II.7) : Principe physique d’un capteur à réseau de Bragg. 

Actuellement l’application des réseaux de Bragg comme capteur a connu un dévelopement 
rapide dans plusieurs domaines parmi lesquels on peut citer : génie-civil, domaine medical, 
les instalations électriques, les structures composites...etc. 


II.5.2.1. Dépendance de la longueur d’onde de Bragg à la température 

Nous avons vu dans la relation (II. 1) que la longueur d’onde de Bragg s’exprime en fonction 
de l’indice de refraction effectif et de la période du réseau. Ces deux grandeurs physiques sont 
aussi dépendantes aux sollicitations externes d’origine thermique, mécanique, hydrostatique, 
électromagnétique et même chimique [7-11]. 

La dépendance de la longueur d’onde de Bragg à la température peut etre déterminée en 
différenciant l’équation (II. 1) comme suit : 


MT) = 2ti(T)a(T) 


(H. 5) 






















Chapitre : II 


Généralités sur les réseaux de Bragg à fibre optique 


La différentiation de l’expression (II.5) par rapport à la température donne : 

Aà b = 2(a^ + 71^W 

b V ST dTJ 


(II-6) 


en divisant l’expression (II.6) par deux coefficients apparaissent : 


Ib 


Ildn 19 a \ 
\n dT ' a dTj 


AT = ( a n + af)AT 


(H-7) 


Où a n — est le coefficient thermo-optique qui correspond à la variation de l’indice de 

réffraction due à une variation de la température AT. 
Pour une fibre optique en silice a n ~ 6.94 * 10 6 /°C. 


a A — ~ 0.5*10 6 /°C. représente le coefficient 

variation de température AT. 


de dilatation thermique du pas due à une 


II.5.2.2. Dépendance de la longueur d’onde de Bragg à la contrainte 

L’application d’une contrainte sur un réseau de Bragg induit une variation de la longueur 
d’onde de Bragg. On peut déterminer ce décalage en différenciant l’équation (II. 1) comme 
suit : 


A4 = 2 (a d f+n d f)Ae (ILS) 

Où As représente la variation de la contrainte appliquée sur le réseau. On considère que le 
réseau de Bragg est soumis à une déformation uni axiale selon Taxe z, et que la matière de la 
fibre optique est homogène et isotrope. En divisant l’expression (II.8) par A, B , on trouve la 
relation suivante : 


Avec : 


AAg _ /I dn 1 9 a \ 

X B \n de a de) z 



(II-9) 


( 11 . 10 ) 
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(“) = —y[Pi2—v(P„ + />i 2 )] (11.11) 

= (1 - T [ p i2 - v ( p ii + p i 2 )])£z (11.12) 

Ab £ 

Où P 12 et P tl sont les coefficients de Pocket et v est le coefficient de Poisson. 


II.5.2.3. Dépendance de la longueur d’onde de Bragg à une pression hydrostatique 

Comme précédemment, on peut déterminer l’influence de la variation de la pression 
hydrostatique sur la longueur d’onde de Bragg en tenant compte que : 

AP 

£ z = -( l-2v)f (11.13) 

Où AP est la variation de la pression hydrostatique et E le module de Young. En substituant 
l’équation (11.13) dans l’équation (11.12), on trouve la relation suivante : 

= -(1 - Y [Pi 2 - v(Pii + Pi 2 )])(l - 2V) f (11.14) 

II.6. APPROCHE THEORIQUE ET MODELISATION DES RESEAUX DE BRAGG 

L’étude et l’analyse de la réponse spectrale du réseau de Bragg est très intéressante. Plusieurs 
approches théoriques ont été utilisées pour modéliser les caractéristiques de ce spectre 
comme la longueur d’onde de Bragg, la réflectivité, la dispersion... etc. Parmi ces méthodes 
théoriques on trouve la méthode de la théorie de Bloch. Elle a été souvent utilisée pour 
décrire le mouvement des électrons dans les semi-conducteurs. Cette méthode a été appliquée 
sur les réseaux de Bragg pour trouver les modes propres de propagation [12]. Une autre 
méthode de modélisation des (FBG), la théorie des modes couplés a été exposée par A.Yariv 
en 1973 [13]. Cette méthode théorique est couramment utilisée, grâce à son taux de réussite. 
C’est un outil approprié pour décrire la propagation des ondes optiques dans un guide d'onde 
avec un indice variant lentement le long de la longueur du guide d'onde et la fibre à réseau de 
Bragg présente ce type de structure [3, 5, 13]. 
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IL 6.1. Théorie des modes couplés 

L’idée principale de la théorie des modes couplés est que le champ électrique du guide 
d’ondes avec perturbation peut être exprimé par une combinaison linéaire des modes de la 
distribution du champ, mais sans perturbation [3, 5, 13]. 

Le champ électrique modale de la fibre optique peut être représenté par : 

E+(x,y,z ) = e+j(x,ÿ)e ±l Pi z j=1.2.3. (11.15) 

Où e ± j(x, y) est l’amplitude du champ électrique transversale du j ieme mode de propagation, 
± représente la direction de propagation et /?, la constante de propagation. En général chaque 
mode de propagation admet sa valeur unique de /? 7 . La propagation de la lumière le long des 
guides d’onde optique repose sur l’usage des équations de Maxwell (chapitre I). 

En notion de la théorie des modes couplés, la composante transverse d’un champ électrique à 
une position z dans une fibre perturbée (fibre à réseau de Bragg) peut être décrite par une 
superposition linéaire des modes guidés d’une fibre non perturbée, et qui est donné comme 
suit : 


Ë t (x, y, z, t) = (x, y, z, t) + E_j(x, y, z, t)) (II. 16) 

En substituant l’équation du champ modal (11.15) dans (11.16), le champ électrique 
E t (x, y, z, t) est écrit par [3,5] : 

Ê t (x,y, z, t) = ü;[riy t '(z)e l ^ z + Aj{z)e~ l Pi z ] e,(x, y)e~ lù)t (11.17) 

Où Aj~ (z) et AJ (z) sont les amplitudes des ondes incidentes et réfléchies respectivement, Pj 
la constante de propagation et êj le champ du mode transverse. Avec E t (x,y, z,t) , est une 
des solutions des équations de Maxwell. 

L’indice de réfraction du réseau de Bragg n(x, y, z) varie selon la position z le long de la fibre 
et s’exprime par : 

2 TC 

n(jx,y,z ) = n(z) = n 0 + ôn 0 + 5n(z)cos[—z + 0(z)] (11.18) 

A 
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Où n 0 est l’indice de réfraction du cœur de la fibre sans perturbation, n — n 0 + ôn 0 la 
valeur moyenne de la modulation de l’indice, en général n 0 » 5n 0 , Sn(z) l’amplitude de 
modulation de l’indice le long de la fibre, 0(z) la phase du réseau et a la période du réseau 
de Bragg. 

La distribution du champ électrique dans le réseau, E t (x,y,z,t ) satisfait l’équation scalaire 
de la propagation d’onde suivante : 

[Vt + k 2 n 2 (x,y, z) — /3 2 ]Ë t (x,y,z,t) — 0 (11.19) 

2n 

Où k = — , représente le vecteur d’onde, /? la constante de propagation dans le vide et A la 

A 

longueur d’onde dans le vide. 

En substituant Le champ électrique E t (x,y,z,t) et l’indice de réfraction n(x, y, z) dans 
l’équation de propagation (11.19) on obtient les équations des modes couplés suivantes : 


^n — i y A+ (lSt , K z \ i(Pm-Pn)z + / V A~ (fZ 1 _ 
1 Lum \ lv mn ' lv mnJ c ' 1 \ lv mn 

— =-iY A + (vt _ JZZ '\ dp m +Pn)Z-iY A~(K t 
d z 1 Aam \ lv mn Iv mnJ^ 1 Lam \ lv mn 


Jf z i(Pm+Pn) z 
lx mnJ c 

-L fA z > \p~^iPm~Pn) z 


( 11 . 20 ) 

( 11 . 21 ) 


Où /é^ n (z) est le coefficient de couplage transversal entre les modes m et n donné par [5] : 

Kmn(z) = |/£ dxdyAs{x, y, z)e mt (x, y)e^ t (x, y ) (11.22) 

As, représente la perturbation de la permittivité, approximativement donné par : 

As s 2 n 0 5n 0 si n 0 » ôn 0 . (11.23) 

Le coefficient longitudinal K est analogue à K^ n , mais en général, /é^ n » K^n pour les 
modes de fibre et ce coefficient-là est donc habituellement négligeable. 


41 



Chapitre : II 


Généralités sur les réseaux de Bragg à fibre optique 


II.6.2. Modélisation des réseaux de Bragg 

L’application de la théorie des modes couplés dans le cas du réseau de Bragg repose sur les 
approximations suivantes [5, 14] : 

■S la fibre est monomode et sans perte, il n’y a pas de modes de propagation en dehors du 
cœur de la fibre, ce qui permet de négliger les modes de gaine, 

■S l’interaction a lieu seulement entre les modes contra-propagatifs dans le cœur, 

■f la différence entre les indices de réfraction du cœur et de la gaine étant très petite. 

La distribution du champ électrique (11.16) selon l’axe de propagation z est décrite, en termes 
de modes fondamentaux incidents et réfléchis, comme suit [5,14]: 

E(x,y,z) — [A + (z)e~ l P z + A~(z)e l ^ z ]e t (x,y) (11.24) 

On remplace l’équation (11.24) dans les équations des modes couplés (11.20) et (11.21) et on 


obtient les équations des modes couplés simples suivantes : 

^- = iâ(z)U(z) + iK(z)V(z) (11.25) 

= -iâ(z)V(z) - iK*(z)U(z) (11.26) 

Où U (z) et V (z) sont respectivement le mode incident et le mode réfléchie, donnés par : 

f/(z) =i 4 + (z)e(‘ 5z- 2 ) (H.27) 

L(z) = 21 -(z)e (_i5z ~ 2 ) (11.28) 


Avec â(z) le coefficient de l’auto-couplage (DC) et K(z) représente le coefficient général 
de couplage (AC) qui permet l’échange d’énergie entre les modes. Ces deux coefficients sont 
très importants pour la simulation de la réponse spectrale d’un réseau de Bragg [5,15]. 

Le coefficient de l’auto-couplage ér(z) est donné par cette relation : 

â = 5 + a-i— (11.29) 

2 dz 
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Où-décrit le « chirp », variation de la période du réseau et 0 la phase. Le désaccord 

2 dz 

« detuning » 8 est représenté comme suit : 


s = p-- = p-p B 

A 

(11.30) 

= 2wn eff (i-A) 

(11.31) 


Où À B — 2n e f jA , est la longueur d’onde conçue pour la réflexion de Bragg. Pour une 
variation minimale de l’indice de réfraction effectif ( 8n e ff ~ 0) on a : 

<T = ySn e// (11.32) 


Le coefficient général de couplage K (z) est donné par cette relation : 

K — K* — jôn(z)g(z)v 

A 


(11.33) 


Où g(z) représente la fonction d’apodisation et v la visibilité de frange. 


II.6.3. Solution des équations des modes couplés 

À partir des équations des modes couplées (11.25) et (11.26) et l’utilisation des conditions au 
limites de la propagation d’un champ optique à travers le réseau, on peut calculer le 
coefficient de réflexion et de transmission d’un réseau de Bragg. 

Considérons donc un réseau de Bragg de longueur L, a son pas, n e ff l’indice de réfraction du 
mode guidé et À B la longueur d’onde de Bragg. La figure (II.8) décrit les conditions aux 
limites de la propagation des modes contra-propagatifs t/(z) et V (z). Le signal incident à 
partir du côté gauche du réseau donne t/(0)=l et L(L) = 0 . À l’aide de ces conditions 
initiales et les équations des modes couplées on obtient le coefficient de réflexion comme 
suit : 
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Figure (II.8) : Réflexion et transmission des modes guidés dans un réseau de Bragg 


Côté gauche : 


{/( 0 ) = 1 
R(0) =? 


(11.34) 


Côté droit : 


f U(L ) =? 

I V{L) = 0 

L’amplitude du coefficient de réflexion r est donnée par : 

_ yfo) 

U( 0) 

La réflectivité du réseau de Bragg R est écrite comme suit : 

R = |r 2 | 


(11.35) 


(11.36) 


(II. 37) 


II.6.3.1. Réseaux de Bragg uniformes 

Pour un réseau de Bragg uniforme, les deux coefficients de couplage K et â sont des 
constantes. Cependant on peut résoudre les équations des modes couplées (11.25) et (11.26) 
analytiquement. En dérivant ces équations et à l’aide des conditions aux limites décrites 
auparavant par les équations (11.34) et (11.35), on obtient: 
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d 2 U 

dz 2 


= (K 2 


d 2 V 

dz 2 


= (K 2 


-â 2 )U 

-â 2 )V 


(11.38) 


Dans ce cas l’amplitude du coefficient de réflexion est obtenue par cette relation [5, 14,15] : 


r = V(0) = 


Ksinh(^K 2 -â 2 L) 


cr sinh (sj K 2 -a 2 L) + isj K 2 - a 2 cosh (yl K 2 -<j 2 L) 


(11.39) 


La réflectivité du réseau R est donnée par : 

sinh 2 (s!k 2 -â 2 L) 


R = r~ = 


(11.40) 


cosh 2 (s!k 2 -â 2 L)~ ° 2 


K 2 


Pour un réseau uniforme, la fonction d’apodisation (g(z) = 1) et la visibilité v prend 
également la valeur 1. À la longueur d’onde de Bragg À B , le coefficient de l’auto-couplage 
(ér=0) et le réseau présente un pic de réflectivité R max suivant cette forme: 

R max =tanh 2 (|^|L) (11.41) 

Les caractéristiques de la réflectivité peuvent être exprimées par les paramètres suivants [16] : 

1. la longueur du réseau L, 

2. la longueur d’onde centrale, À B , 

3. l’indice effectif n e jf, 

4. l’amplitude de modification de l’indice effectif, An e ^, 

5. le pic de réflectivité R max . 

La réflectivité d’un réseau peut atteindre 1 si on augmente la longueur du réseau ou la 
modulation de l’indice de réfraction An. Les figures (II.9.a), (II.9.b), (II.9.c) et (II.9.d) 
illustrent l’allure de la réflectivité d’un réseau unifonne simulée pour différentes valeurs de la 
longueur du réseau. Les valeurs de simulation sont : L= (5mm, 7mm, 10mm, 15mm), 
A s =1550mm, n e jf =1.47, An e ^=l e ' 4 . 
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Figure (II.9) : Réflectivité d’un réseau uniforme pour différentes valeurs de L [16]. 


À partir de la figure (II. 9) on constate qu’avec l’augmentation de la longueur du réseau on 
obtient une réflectivité plus importante, mais on observe aussi que le nombre des lobes 
secondaires augmente aussi. Ces lobes peuvent être éliminés par l’utilisation des réseaux de 
Bragg apodisés. 


ûrt=0 0005 




Figure (11.10): Réflectivité d’un réseau uniforme pour différentes valeurs de An [17]. 
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La figure (II. 10) présente l’allure de la réflectivité d’un réseau de Bragg unifonne ayant les 
paramètres suivants L=4mm, À B = 1309.45 nm, n e ff =1.441, simulée en fonction de la longueur 
d’onde pour différentes valeurs de An e jf (0.00005, 0.0001, 0.0003, 0.0005). On remarque que 
nous obtenons une réflectivité plus importante et une largeur spectrale plus grande en 
augmentant la valeur de A n eJ f . Cette conclusion assure la possibilité de concevoir des réseaux 
de Bragg courts qui fournissent quand même de très bonnes réflectivités en diminuant À n e g . 


II.6.3.2. Réseaux de Bragg non uniformes 

La résolution des équations des modes couplés (11.25) et (11.26) analytiquement pour un 
réseau de Bragg non uniforme est très difficile. Plusieurs méthodes numériques ont été 
employées [18,19]. Mais la plus commune et la plus utile pour simuler les spectres de 
réflexion ou transmission d’un réseau non uniforme est la méthode des matrices de transfert 
(TMM, Transfer Matrix Méthode). Cette méthode numérique a été décrite spécifiquement 
pour la première fois par Yamada [20]. 

La méthode des matrices de transfert consiste à diviser le réseau non uniforme de longueur L 
en M segments uniformes de longueur L î; figure (11.11). A chaque sous réseaux i correspond 
une matrice de transfert F,. Les termes de cette matrice d’ordre 2 sont décrits comme 
suit [5,16]: 



cr 


-i — smh(y B dz) 

Yb 


cosh (y B dz) - i —sinh (y B dz) 

Yb 

i — sinh (y B dz) cosh (y B dz) + i — sinh(y B dz) 

Yb Yb 


(11.42) 


Où â est le coefficient de l’auto-couplage donné par la relation (11.29), K le coefficient 
général de couplage défini par la relation (11.32) et y B est décrit comme suit : 


y B =^K 2 -a 2 (K 2 >â 2 ) 

(11.43) 

y B = isjà 2 - K 2 (K 2 < â 2 ) 

(11.44) 
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Figure (11.11) : Principe de la méthode des matrices de transfert. 


On peut ensuite calculer le produit des matrices individuelles pour trouver la matrice totale 
présentée par cette relation : 


Où 


~U(L) 

= \F lot ] 

'£/( 0)" 

V(L) 


_V(0)_ 


.[i 7 ;]x...[F 2 ]x[F 1 ] 



(11.45) 


(11.46) 


Avec l’utilisation des conditions aux limites qui sont décrites par les équations (11.34) et 
(11.35), on obtient enfin le coefficient de réflexion présenté par : 


r = --41 (H.47) 

F 

1 22 

Et le coefficient de transmission donné par : 

t = ^r (H.48) 

Ml 


II. 6.4. Calcul du délai de groupe et de la dispersion 

La dispersion de vitesse de groupe induit dans un guide d’onde optique perturbe la 
propagation du signal dans les réseaux de télécommunication. Ce phénomène se traduit par 
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l’élargissement temporel des impulsions lumineuses constituant le signal, ce qui a pour effet 
d’introduire des distorsions dans le signal transmis. Il est donc intéressant de caractériser les 
composants optiques en termes de dispersion. On peut calculer le délai de groupe et la 
dispersion d’un réseau de Bragg à partir de la phase du coefficient de réflexion [15]. 


Le délai de groupe r r (ps) est défini par cette relation : 

_ d6 r A 2 d6 r 

r dœ 2nc dA 


(11.49) 


Ou 6 r est la phase du coefficient de réflexion (r). 

La dispersion d r (en ps/nm) est donnée par cette relation : 

, _ dx r _ 2r r A 2 d 2 d r _ 2nc d 2 0 r 
r dA A 2nc dA 2 A 2 dco 2 


(11.50) 


II.7. CONCLUSION 

Dans ce chapitre nous avons présenté la structure d’un réseau de Bragg photo inscrit dans une 
fibre optique ainsi que les différentes techniques d’inscription utilisées. Ensuite nous avons 
décrit les deux paramètres caractéristiques du réseau de Bragg, le pas et l’indice de réfraction 
ainsi que leur dépendance aux sollicitations externes d’origine thermique, mécanique et 
hydrostatique. Cette sensibibité représente un atout dans le domaine des capteurs à réseaux 
de Bragg. Les applications des réseaux de Bragg dans le domaines des telecomunications ont 
été également présentés. 

Une étude théorique sur la propagation des ondes lumineuses dans la structure FBG a été 
décrite par les équations des modes couplés. Ces équations peuvent etre résolu 
annalytiquement pour un réseau uniforme, cependant pour un réseau non uniforme, 
l’utilisation d’une méthode numérique s’avère nécéssaire . Dans ce chapitre nous avons décrit 
la méthode des matrices de transfert (TMM) pour simuler et analyser le spectre de 
réflectivité d’un réseau non unifonne et cette intéressante approche numérique a démontré son 
efficacité. 
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III.1. INTRODUCTION 

L’étude de la propagation des ondes de types solitons dans les guides d’ondes optiques constitue 
un sujet de recherche fondamental et très intéressant. Ce phénomène est apparu dans les 
dernières années du vingtième siècle comme une technique extrêmement prometteuse pour 
transmettre sur des fibres optiques des débits très élevés sur de très grandes distances. 

Un soliton optique est une impulsion lumineuse pouvant se propager sans déformation sur des 
distances extrêmement grandes (théoriquement infinies). Il reste stable au cours de sa 
propagation et ceci grâce à la compensation mutuelle de la dispersion chromatique et l’effet non 
linéaire (Kerr) [1]. La dispersion a pour effet l’étalement de l’impulsion par contre l’effet Kerr 
entraine l’effet inverse. 

Comme de nombreux guides d’ondes optiques, les réseaux de Bragg fibrés sont susceptibles de 
donner lieu à la propagation d’ondes solitaires. Ce phénomène a été mentionnée pour la première 
fois par Chen et Mills en 1987 [2], avec la propagation d’un soliton particulier appelé « soliton 
de Bragg ». De même que le soliton conventionnel dans la fibre optique, le soliton de Bragg 
résulte des deux effets combinés, l’effet Kerr (auto-modulation de phase) et la dispersion de 
vitesse de groupe induite par le réseau. Ce soliton peut se propager sans déformation dans le 
réseau avec une fréquence située dans la bande interdite et une vitesse pouvant avoisiner celle 
de la lumière. 

Nous abordons ce chapitre par un rappel historique sur le concept soliton. Les différents types 
de solitons optiques pouvant se propager dans une fibre optique seront aussi présentés. Une 
partie du chapitre sera consacré à l’étude théorique de la propagation des solitons de Bragg dans 
la structure FBG. Aussi, l’utilisation de la méthode numérique SSFM dans notre contribution 
concernera l’étude de la dynamique non linéaire des solitons optiques. 


III.2. RAPPEL HISTORIQUE 

Le soliton est une onde solitaire qui se propage dans les milieux non-linéaires et dispersifs. 
Cette onde spécifique, possède des propriétés de stabilité remarquables, (localisation spatiale, 
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conservation de la vitesse de propagation et de l’énergie). La première observation d’un soliton 
fut réalisée par l’ingénieur écossais John Scott Russell en 1834 [3]. Alors qu’il se promenait à 
cheval le long d’un canal dans la région d’Edimbourg, un bateau qui était tiré le long de ce canal 
s’arrêta brusquement près de lui. Ce bateau mit en mouvement dans le canal une masse d’eau 
que J.S. Russell tenta de suivre. Il observa alors que cette masse d’eau prit la fonne d’un paquet 
d’eau rond qui se propageait dans le canal sans changement de fonne ni variation de vitesse. Il 
perdit ensuite la trace de cette « grande onde solitaire » après l’avoir suivi pendant un ou deux 
miles, figure (III. 1). 



Figure (III.l) : Simulation de l’observation de J.S. Russell [3]. 


John Scott Russell a été tellement fasciné par ce phénomène étrange, qu’il a consacré une dizaine 
d’années de sa vie à étudier le concept de cette vague solitaire. Cependant il fallut attendre 1895 
pour que l’équipe de Kertweg et Vries décrit ce phénomène théoriquement par le développement 
d’une équation non-linéaire dite de « Kertweg et de Vries » (KDV). La forme simple de cette 
équation est donné par : 

ÔA{Z,t) i 6A Ô A(Z,t) l ^*> -0 (III.l) 

dt 8z dz 3 

Où A(z, t) représente l’amplitude du champ hydrodynamique qui se propage à la surface de l’eau 
selon l’axe z. 
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Cette équation admet comme solutions analytiques des solitons dont celles de Russell et ceux 
sont Zabusky et Cruskal les premiers a les avoir trouvé en 1965 par l’utilisât ion de la méthode 
transformation inverse [1,4]. Ces solutions présentent une allure de type séquente hyperbolique 
comme suit: 


A(z,t ) = — sec h' 
2 


— (Z - vr) 
2 


(III.2) 


Où v est la vitesse du soliton, et sech est la fonction séquente hyperbolique. 


D’un point de vue fondamental, la formation des solitons dans les milieux non linéaires et 
dispersifs est dû à un exact équilibre entre la dispersion (ou la diffraction) qui tend à l’expansion 
de l’onde localisée, et l’effet non-linéaire qui contre balance l’effet dispersif. Cette 
caractéristique est unique : cela implique qu’à travers l’effet non-linéaire, l’onde induit 
simultanément un puit de potentiel et se capture elle-même dans son propre potentiel induit. 
Encore plus fascinant et en se propageant, les solitons peuvent interagir avec d’autres solitons et 
possèdent donc à ce titre les mêmes propriétés qu’une particule. Cependant les solitons ou les 
ondes solitaires se manifestent donc naturellement dans la plupart des systèmes non-linéaires et 
sont même devenus indispensables pour la description de nombreux phénomènes physiques tels 
que [1-7]: 


■ la propagation de certaines ondes en mécanique des fluides, 

■ les ondes localisées dans les plasmas astrophysiques, 

■ la propagation de faisceaux lasers dans les fibres optiques, 

■ les phénomènes de transport de charge dans les polymères conducteurs, 

■ les modes localisés dans les cristaux magnétiques, 

■ les dynamiques de macromolécules biologiques comme l’ADN et les protéines. 


La propagation des ondes solitaires dans les fibres optiques est couramment modélisée par 
l’équation non linéaire de Schrôdinger. Cette équation peut décrire tous les effets linéaires et non 
linéaires qui affectent l'impulsion au cours de sa propagation dans la fibre optique (voir chapitre 
I). La forme simple de cette équation est donnée par : 
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dA . B , d 2 A 

-h l - — 

ôz. 2 ÔT 2 


iy A A = 0 


(III.3) 


L’équation (III.3) inclut le terme de dispersion de deuxième ordre B 2 , et la non linéarité de type 
Kerr y. Cette équation représente un modèle complètement intégrable, où il est possible d’obtenir 
des solutions analytiques de types soliton. En 1972, Zakharov et Shabbat, ont été les premiers à 
résoudre cette équation et ceci grâce à la méthode de la diffusion inverse [8], Ils ont trouvés 
deux types de soliton : pour Po>0, (dispersion nonnale), la solution prend la forme d’une 
tangente hyperbolique qui représente un soliton sombre. Mais dans le cas où [L<0 (dispersion 
anormale), la solution prend la forme d’une séquente hyperbolique qui représente un soliton 
brillant. En 1973, Hasegawa et Tappert ont démontré numériquement la stabilité des impulsions 
pico-seconde de type soliton dans les fibres optiques [9,10]. Ensuite, en 1980, le soliton optique 
fut démontré expérimentalement par Mollenauer et ses collaborateurs [11]. De nos jours, de 
nombreux laboratoires de recherche en transmission de l’infonnation travaillent sur la 
dynamique et la stabilité des solitons dans les fibres optiques et les structures périodique tels que 
les réseaux de Bragg. En effet, le soliton ayant la capacité de se propager en gardant ses 
propriétés, il peut permettre de concevoir des systèmes optiques plus performants et donc 
d’améliorer grandement les technologies de l’information et de la communication. 


III.3. DIFFERENTS TYPES DE SOLITONS OPTIQUES 

En effet, toute impulsion ou paquet d’ondes a une tendance naturelle à s’étaler durant sa 
propagation dans un milieu. En optique, une onde localisée dans l’espace ou dans le temps peut 
subir un étalement, soit de son enveloppe temporelle ou bien de ses dimensions spatiales ou 
même des deux simultanément. 

Pour une impulsion temporelle, l’étalement est dû à la dispersion chromatique, par contre 
l’étalement d’une impulsion spatiale est dû à la diffraction. En régime non linéaire, les effets non 
linéaires jouent un rôle principal dans la compensation de ces étalements et donnent lieu ainsi à 
l’apparition des solitons. Dans ce concept, il en résulte trois types de soliton optique : les solitons 
temporels, les solitons spatiaux et les solitons spatiaux temporels ou bille de lumière. 
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III.3.1. Solitons temporels 

La balance entre l’effet de dispersion chromatique et l’auto-modulation de phase induit par 
l’effet Kerr, est la clé de la formation des solitons temporels. En effet, la dispersion chromatique 
induit un élargissement temporel de l’impulsion, en opposition, l’auto-modulation de phase 
(SPM) induit l’effet inverse (compression de l’impulsion) [9,10]. 

L’auto-modulation de phase qui dépend de l’indice non linéaire « 2 , et de l’intensité 7, se traduit 

par un déphasage non linéaire A</> SPM = n 2 k 0 IL de l’impulsion après une propagation sur une 

distance L. Ce déphasage induit un décalage de fréquence _ -\ œ . Il est donc possible 

dt 

d’obtenir une impulsion soliton si le décalage de fréquence compense exactement la dispersion 
chromatique. 


III.3.1.1. Soliton brillant 

En régime de dispersion anormale ([L<0), les hautes fréquences sont plus rapides que les basses 
fréquences. En opposition l’effet non linéaire (SPM) va produire un décalage de fréquence se 
traduisant par le ralentissement des fréquences élevées et l’accélération des fréquences basses. 
Ainsi, sous certaines conditions, un équilibre peut s’instaurer et l’effet Kerr maintient alors en 
phase les composantes fréquentielles qui se seraient étalées par dispersion. 


L’impulsion qui en résulte, est appelée soliton brillant: elle est couramment rencontrée dans les 
fibres optique qui possèdent un régime de dispersion anormale. Théoriquement, elle correspond 
à une famille de solutions analytiques de l’équation de Schrôdinger (III.3), avec un profil en 
séquente hyperbolique et une phase plane. La solution d’un soliton fondamental brillant se 
présente comme suit [11]: 


u(Z,t ) = sec h(r) exp (/Z / 2) 


(III.4) 


Où Z et x, sont la distance et le temps normalisés respectivement. 


L’amplitude normalisée u est définie comme suit: 

u = N^L = J^A 
V ° 0 


(HL 5) 
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Avec P Q , la puissance crête de l’impulsion et N l’ordre du soliton donné par : 

A^ 2 = —= (III. 6) 

Lnl |-®21 

Le soliton fondamental (A=l), se propage sans déformation par rapport aux soliton d’ordre 
supérieur. De ce fait, il est le plus attractif pour les transmissions optiques. La figure (111.2) 
présente le comportement d’un soliton fondamental dans une fibre optique. 



Figure (III.2) : Evolution d’un soliton fondamental brillant dans une fibre optique. 


III.3.1.2. Soliton noir 

En régime de dispersion normale (P2>0), les fréquences élevées sont cette fois plus lentes que les 
fréquences basses. Pour compenser l’effet d’étalement du à la dispersion chromatique, il faut 
inverser le signe de la variation temporelle de la loi d’auto-modulation de phase. Ainsi, 
l’équilibre peut être instauré et l’impulsion qui en résulte dans ce cas est appelée « soliton noir ». 
C’est une impulsion particulière qui possède deux fronts de phase plane de part et d’autre du 
creux, d’un profd de forme tangente hyperbolique. L’équation (III.7) présente le profil d’un 
soliton fondamental noir [11]: 


u(Z,t) = tanh(r)exp(/Z) 


(III.7) 
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La figure (III.3) présente l’évolution d’un soliton fondamental noir. 



Figure (III.3) : Evolution d’un soliton fondamental noir dans une fibre optique. 


III.3.2. Solitons spatiaux 

Les solitons spatiaux sont des faisceaux optiques dont la diffraction naturelle a été compensée 
par l’effet non linéaire du milieu de propagation sensible à l’intensité. La diffraction induit un 
élargissement du faisceau, par contre l’effet non linéaire agit comme une lentille convergente ou 
divergente selon que l’indice est augmenté ou diminué. Cet effet va permettre l’auto focalisation 
du faisceau durant sa propagation. Tout comme dans le domaine temporel, quand les deux effets 
l’auto-focalisation et la diffraction naturelle se compensent exactement, l’observation d’un 
soliton spatial sera possible. Ces solitons présentent eux aussi un intérêt pour les 
télécommunications [13]. 


III.3.3. Billes de lumière 

Au carrefour des solitons spatiaux et temporels, les solitons spatio-temporels viennent compléter 
la gamme déjà riche en diversité des solitons. La propagation d’une impulsion optique intense et 
focalisée est gouvernée par l’entrejeu entre la diffraction, la dispersion de vitesse de groupe et les 
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différents processus non linéaires présents. Sous certaines conditions, les non linéarités peuvent 
être utilisées pour compenser à la fois la diffraction et la dispersion chromatique, et ainsi 
produire simultanément un soliton dans l’espace et dans le temps. Des études expérimentales 
récentes ont montré que les matériaux possédant une forte non-linéarité avec un indice de 
réfraction non linéaire d’ordre quatre négatif devaient permettre cette propagation solitonique 

[14.15] . Au-delà de leur intérêt scientifique, ces véritables « boules de lumière » constituent une 
source d’applications potentielles notamment en traitement tout optique de l’information 
(commutation ultrarapides, contrôle de trajectoire, opération logiques) à trois dimensions 

[14.16] , 

III.4. SOLITONS BRAGG 

Les solitons Bragg sont réalisés par l’équilibre de la balance entre la dispersion du réseau de 
Bragg et la non-linéarité de type Kerr de la fibre. En effet la périodicité de l’indice de réfraction 
induit une bande photonique interdite (PBG), dont la position de cette bande sur l’axe fréquentiel 
dépend de l’indice de réfraction du milieu. Dans la présence d’effet Kerr non linéaire, l’intensité 
du champ optique est suffisante pour modifier l’indice de réfraction. Cependant une fréquence 
située dans la bande interdite à faible puissance ne le sera pas nécessairement à plus forte 
puissance. Ainsi, une impulsion intense et dont la fréquence se trouve dans la bande interdite 
pourra être confinée dans le réseau. Dans les ailes de cette impulsion, l’intensité est faible et il y 
aura réflexion de Bragg, par contre, en son centre, les ondes régressives et progressives pourront 
se propager, empêchant l’impulsion de disperser dans la structure, et permettre ainsi la formation 
d’impulsion solitonique (Gap soliton), qui peuvent se propager sans déformation. Avec cette 
configuration, le soliton de Bragg est formé par la superposition de deux ondes de fréquences 
différentes, l’onde progressive et l’autre régressive. Si les intensités des deux sont égales, 
l’enveloppe de l’impulsion restera immobile dans le réseau, c’est le cas d’un soliton de Bragg 
stationnaire [2]. Dans le cas où les intensités sont différentes, l’enveloppe de l’impulsion se 
propagera dans la direction de la plus intense des deux, avec une vitesse comprise entre zéro et la 
vitesse de groupe de ces ondes en l’absence de réseau [17]. 
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En dehors de la bande interdite, les solitons de Bragg peuvent aussi exister, mais dans ce cas, 
leur formations s’explique par un équilibre entre la dispersion du réseau et la non-linéarité, en 
parfaite analogie avec les solitons temporels dans les fibres optiques 

L’existence des solitons dans les réseaux de Bragg a été mentionnée pour la première fois en 
1987 par Chen et Mills [2], Par simulation numérique, ils ont démontré la possibilité d’auto¬ 
localisation de la lumière dans le réseau de Bragg sous forme d’un soliton particulier dit « soliton 
de Bragg ». Ensuite, en 1988, Sipe et Winful puis Destterke et Sipe [17,18], ont montré que, dans 
une certaine limite, la propagation des impulsions optiques dans les FBG pouvait être modélisée 
par un système d’équations aux modes couplés linéaires en régime linéaire, et par un système 
d’équations aux modes couplés non linéaires en régime non linéaire. Bien que non intégrable, ce 
système d’équations non linéaires possède des solutions analytiques localisées. En 1989, Aceves 
et Wabnitz [19], Christodoulides et Joseph [20], sont les premiers à avoir réussis l’obtention de 
ces solutions analytiques. Peux après en 1993, Feng et Kenbuhl [21], ont trouvé des solutions 
analytiques de type soliton brillant ou noir, mais dont l’amplitude du champ ne s’annule pas à 
l’infini. 

Expérimentalement, plusieurs travaux ont été réalisés pour démontrer et observer la formation 
des solitons de Bragg. En 1998 Taverner et ses collaborateurs [22] ont rapporté la première 
observation de la génération du soliton Bragg dans la zone de la bande interdite. Ensuite, en 
1999 Eggleton et al [23], ont réalisé avec succès la première observation du soliton Bragg en 
dehors de la bande interdite où le réseau est transmissive mais fortement dispersif. 

Récemment, les applications des solitons Bragg dans la réalisation des systèmes bistables ont été 
extrêmement étendues. On les retrouve dans la compression des impulsions, la limitation, les 
opérations de logique et finalement, vu que la vitesse du soliton Bragg peut être arbitrairement 
faible, il pourrait être utilisé comme mémoire temporaire. 


III.4.1. Equation des modes couplés non linéaires 

La compréhension théorique de la propagation non linéaire des impulsions optiques dans les 
réseaux de Bragg fibrés repose sur la théorie des modes couplés. Cette dernière est largement 
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appliquée pour décrire la propagation linéaire dans les structures périodiques telles que les 
réseaux de Bragg (voir chapitre II). Mais elle présente la simplicité et la facilité d’utilisation en 
présence des non linéarités. 


Afin d’établir le système d’équations des modes couplés non linéaires, commençons par 
l’équation d’onde qui décrit la propagation de la lumière dans les fibres optiques (voir chapitre 
I), donnée par cette relation [24 ,25]: 


Où la polarisation P induite dans la fibre optique peut être exprimée par cette relation : 


(HL 8) 


P = ^oZ l E + £ 0 s NL É 


'OA 1 
3 31*712 


£nl=~X j E 


(III.9) 
(III. 10) 


En présence d’un réseau de Bragg uniforme le champ électrique prend la forme suivante : 

É(r, t) = | F(x, y)(A + (z, t)e~ i(û¥ - kl>z) + A (z, t)e~ l<û¥+k » z) ) ( m - 1 1) 

Où A + et A sont les enveloppes des deux ondes progressive et régressive respectivement. 

Avec : 


n 


k B = — et (o Q qui représentent respectivement le vecteur d’onde et la fréquence centrale. 
A 

Dans le domaine fréquentiel, le champ électrique satisfait l’équation de Helmholtz suivante : 


V 2 Ë + £(co)^Ê = 0 


(III. 12) 


Avec : 


Ë(r, co) = F(x, y)(Â + (z, co)e i(k ° z) + A (z, co)e- i(k ° z) ) 

Où £(co) est la fonction diélectrique tel que : 


(III. 13) 


£(z) = l + X il \z) + £ 


NL 


(III. 14) 


En utilisant la relation entre la fonction diélectrique et l’indice de réfraction 
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s(a>) = n 2 (co) (III. 15) 

Où la modulation de l’indice de réfraction en régime non linéaire est donnée par : 

n(co,z) = n(co) + Sn + n 2 \E\ (III. 16) 

On introduit les équations (III. 13) et (III. 16) dans l’équation d’onde (III. 12), et avec quelques 
simplifications et quelques transformations algébriques on obtient dans le domaine fréquentiel 
les équations non linéaires des modes couplés suivantes [23, 24] : 


ÔA ~ ~ 

—— = i[S(a>) + A B]A + + ixA_ 

dz 

(III. 17) 

ÔA ~ ~ 

-— = i[S(a>) + A B]Â + ixA + 

dz 


Où ù(eu), est le désaccord de fréquences défini par : 

S(ûO = ( n /c)(<0 - C o B ) = 5(m) - k B (III. 18) 

Et k représente le coefficient de couplage linéaire, qui pennet l’échange d’énergie entre l’onde 
progressive et l’onde régressive. Pour un réseau de Bragg uniforme, k est donné par cette 
relation : 


2nSn 



(III. 19) 


Le terme AB où AB = B (m) — B(o)) : représente la contribution non linéaire à la propagation. 
S (eu) , est la constante de propagation, développée en série de Taylor autour de la pulsation 
centrale de l’impulsion co 0 com 

p{(o) = po + A («J-(O 0 ) + ^/3 2 (oj-(O 0 Y +^-/3 3 {co-co o y +A-/3 4 (a>-a> 0 ) 4 +..(111.20) 

2 6 24 
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Le système d’équations des modes couplés (III. 17) peut se convertir dans le domaine temporel. 
En remplaçant (eu — co 0 ) par l’opérateur différentiel (ij-) , le système d’équations des modes 
couplés non linéaires s’écrit comme suit: 


ÔA + dA + iB , ô A 
+ +B ,—- + - 


& 

ÔA 


dt 


2 dt 


2 + ^ A + = iôA + + ihA + z'L(|A + 1 2 + 2| A_ | 2 )A + 


(III.21) 


+ B , 


8A iB , ô 2 A 


a 


„ H—A + ïkA, + iT( A + 2A. )A 

5r 2 ' + 1 1 +l “ 


ôz 1 dt 2 

Les paramètres /?,,/? 2 représentent respectivement l’inverse de la vitesse de groupe v g et la 
dispersion de la vitesse de groupe (GVD), a le coefficient de perte et r le coefficient de 
couplage non linéaire donné par cette relation : 


r = 


n 2 ù > 0 


cA 


eff 


(III.22) 


On voit que la non-linéarité fait apparaitre deux termes dans chacune des équations du système. 
Les termes non linéaires précédés de 2 sont la modulation de phase croisée (XPM), et les autres 
sont dites auto-modulation de phase (SPM). Les deux équations (III.21) sont les équations des 
modes couplés non linéaires représentant la base sur laquelle la théorie du soliton de Bragg a été 
élaborée. 


III.4.2. Solutions linéaires 

La propagation linéaire dans le réseau de Bragg fibré est modélisée par les équations des modes 
couplés linéaires. Ces équations ont été obtenus en négligeant la contribution des effets non 
linéaires dans le système (III. 17), c’est à dire (A B = 0) [24] : 


' dÂ + 
dz 


iô{co)A + + ihA_ 


cÀ 

~dz 


iô{co)A_ + ihA + 


(III.23) 
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En effet l’équation des modes couplés (III.23) peut être résolue facilement dans le domaine de 
fourrier. La solution générale de ces équations prend la forme suivante : 

A + = Æj exp (ik(co)z) + a 2 exp (-ik(co)z) (III.24) 

A = b x exp (ik(co)z) + b 2 exp (-ik(co)z) (III.25) 

Où les coefficients constants ai, a 2 , bi, bo, sont déterminés par les conditions initiales. En 
injectant les solutions (III.24) et (III.25) dans le système d’équation (III.23), on trouve que les 
constantes ne sont pas toutes indépendantes et satisfont les quatre relations suivantes : 

(, k - ô)a x = td> x {k + ô)a 2 = —td> 2 (III. 26) 

(k — ô)b 2 = m 2 (k + ô)b x = -m 1 

Pour que le système d’équation (III.26) possède des solutions non triviales, il faut 
nécessairement que k vérifie la relation de dispersion linéaire suivante : 

k = ±V S 2 — k~ Q = (oo — co B ) = + — -J k~ +k 2 (III.27) 

n 

La variation de la courbe de dispersion Q(k) est illustrée sur la figure (III.4). 



Figure (III.4) : Courbe de dispersion, montrant une bande interdite PBG. 

À partir de la relation de dispersion, on constate que le couplage entre les modes de propagation 
dans le réseau de Bragg ouvre une bande interdite PBG centré à fl — 0 . La largeur de cette 

c 

bande est égale à A fl — 2 k - . 

° n 


V 
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III.4.3. Solutions non linéaires (Gap soliton) 

Lorsque la réponse non linéaire de la fibre est prise en compte, la propagation des ondes dans les 
FBG est profondément modifiée. Des solutions solitons peuvent exister à des fréquences proches 
mais en dehors de la bande interdite. Bien plus, des structures d’ondes solitaires peuvent encore 
se propager dans la bande interdite dite Gap Soliton. 

En 1988 Sipe et Destterke [18] ont montré que la propagation de ces solitons est décrite par une 
équation de Schrôdinger non linéaire (NLSE) standard, qu’on peut dériver à partir des équations 
aux modes couplés (III.20) comprenant des termes représentant la contribution non linéaire à la 
propagation (XPM, SPM). Ce système d’équations non linéaires possède des solutions 
analytiques localisées au sein de la bande interdite du réseau. Aceves et Wabnitz [19], 
Christodoulides et Joseph [20], sont les premiers à avoir réussi l’obtention de ces solutions 
analytiques. 

Le point de départ de leur raisonnement consiste à rechercher des solutions du système (III. 20) 
au moyen d’une transformation appropriée des solitons du modèle de Thirring massif [25]. Ces 
équations se réduisent au modèle de Thirring massif en l’absence des termes d’auto-modulation 
de phase [25]. Les solutions obtenues sont des ondes solitaires qui se propagent sans 
déformation. En revanche, les simulations numériques montrent que le comportement de ces 
solutions est propre à celui des solitons et est généralement utilisé comme critère pour les définir 
[26], 

Ces ondes solitaires ou (Gap Solion) prennent la forme suivante [19, 23, 24] : 


Avec : 


A + (z, t) = a + sec h(J; - i 
A_{z,t ) = a sec/7(<^ + 


= ± 


1 ± v 
1 + v 


k(\-v 

r(3-v 2 ) 


îsin^ 


(III.28) 


(III.29) 


z-V G t . 

, /rsnw 

(IH.30) 


66 





Chapitre : III 


Théorie des solitons. 


6 = — k cos y/ - ——r tan 1 fjcot(u/72)|coth(£)]. 

3-v 2 U V ' A 1 (III.31) 

Le paramètre v des solutions (III.29) peut prendre toutes les valeurs du rang : |v | <1 et détermine 
la vitesse d’ensemble du soliton. Le paramètre V est aussi un paramètre libre pouvant prendre 
toutes les valeurs comprises dans la plage :0 < V < n .11 détermine l’amplitude et la largeur 
du soliton. 


III.4.4. Relation entre N CME et NLSE 


Nous avons vu précédemment que la propagation non linéaire des impulsions optiques dans un 
réseau de Bragg fibré est modélisée par les équations des modes couplés non linéaires (111.21). 
Dans certaines conditions, ces équations peuvent se réduire par des méthodes mathématiques à 
l’équation non linéaire de Schrôdinger (NLSE) (III.3) [24]. La première condition est que 
l’intensité de l’impulsion d’entrée P 0 est plus faible, assez que n 2 \E\ 2 « 8n . Cette condition est 

1 

équivalente à yP 0 « K ou kL nl » 1 avecL NL , la longueur non linéaire donné par L nl — — . 

yPo 

La deuxième condition est que l’effet de dispersion d’ordre trois B 3 induit par le réseau est 
négligeable. Dans ce cas la propagation des impulsions dans le réseau est régie par l’équation 
non linéaire de Schrôdinger suivante [24] : 


dA . B 2 8 ô 2 A 
ôz 2 dT 2 


iy g A~A = 0 


(III.32) 


Où A est l’amplitude de l’enveloppe de l’onde progressive A + et de l’onde régressive A_, 
B 2 ,y g sont respectivement la dispersion de la vitesse de groupe et le coefficient non linéaire du 
réseau, donnés par les relations suivantes: 

q-v’)* < IIU3 > 

sgn (f)v g 2 icv 3 

r g =(^-)r (in. 34) 

* 2v 
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Nous constatons que les deux paramètres B®,y g ne sont pas des constantes mais dépendent du 
paramètre v. Dans la bande interdite où v = 0, l’équation (III.32) n’est pas valable et n’admet 
pas de solutions. Cependant, près mais endors de celle-ci,(v -> 1) , l’équation (III.32) devient 
intégrable et peut être résolue facilement par la méthode de la diffusion inverse . On assiste à 
l’existence de solitons fondamentaux ou d’ordres supérieurs dans le réseau comme dans la fibre 
et qui sont de la forme suivante [24] : 


u(Z, t) = ?7 sec h[rj(r - r s + £Z)]exp [/(// 2 - s 2 )z / 2 - isz + i(j) s ] 


(III. 35) 


Où Z = z/L d et r — T/T 0 sont la distance normalisé et le temps normalisé respectivement, 
u = yjLoYgA est l’amplitude normalisée, L D = 7’ 0 2 /|5^| est la longueur de dispersion. 
rj, e, t s et 0 S '- sont des paramètres arbitraires qui représentent l’amplitude, la fréquence, la 
position et la phase du soliton respectivement. 


III.5. CARACTERESATION DES SOLITONS OPTIQUES 


Afin d’étudier numériquement l’influence des différents effets qui affectent la propagation des 
solitons dans une fibre optique monomode, on utilise la méthode Split-Step-Fourier (SSFM), du 
fait de sa vitesse d’acquisition et de sa précision dans les calculs [11]. Cette méthode numérique 
a été présentée dans le Chapitre I, (Pag 24). 


La propagation des solitons dans la fibre optique est modélisée par l’équation non linéaire de 
Schrodinger suivante [11]: 


SA a A B. d 2 A B , d 3 A 

-1 —A +1 --- 1 -T 

dz 2 2 ÔT 2 6 dT 3 


iy A A = 0 


(III.36) 


Cette équation inclus l’effet de dispersion de deuxième ordre B 2 , l’effet de dispersion de 
troisième ordre fi? l’atténuation linéique a et l’effet non linéaire y. 
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III.5.1. Influence des pertes 

Au cours de sa propagation dans la fibre optique, l’impulsion lumineuse subit des pertes se 
traduisant par une diminution de l’intensité de l’onde. Ces pertes définies par le coefficient 
d’atténuation linéaire a (dB/km) sont majoritairement induit par l’absorption et par la diffusion 
Rayleigh. 

Pour étudier l’impact des pertes sur la propagation d’un soliton fondamental dans une fibre 
optique monomode on choisit les paramètres suivants [11.27]: 

B 2 = -20 (ps 2 /km), B3=0.1 (ps 3 /km), a=0.14 (dB/km), y = 2.43 (1/kmW), To=20ps, k=1.55pm. 
P 0 = 0.02 W. 

Où l’amplitude de l’impulsion d’entrée est donnée sous la forme : 

u (0, x)=sech (z). t — — . (ffl.37) 

T o 

Les résultats de simulation sont présentés par la figure (III.5). 



Figure (III.5) : Evolution d’un soliton fondamental dans une fibre optique, (a) sans perte, 

(b) avec perte. 


7 ^ 
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À partir de la figure (III.5.a) qui représente la propagation d’un soliton fondamental dans une 
fibre optique sans perte c’est à dire (a=0 dB/km), on constate que l’amplitude du soliton n’a pas 
changé et le soliton reste stable et se propage sans distorsion, ceci résulte de la compensation 
mutuelle entre l’effet de dispersion chromatique et l’effet non linéaire de Kerr. 

La figure (III.6.b) montre le cas où la fibre présente une perte (a=0.14 dB/km), on constate que 
l’amplitude du soliton diminue le long de la fibre. Ceci signifie que les pertes sont un facteur 
limitatif qui réduit la puissance du soliton et par conséquent limite les performances de la 
transmission optique par soliton sur les longues distances. 

Pour résoudre ce problème, on installe périodiquement le long de la fibre des amplificateurs à 
fibre dopées erbium (EDFA) [11] qui ont montrées à ce jour leurs efficacités surtout dans le 
domaine des télécommunications par fibre optiques. 


III.5.2. Interaction entre les solitons 

Le paramètre le plus important dans la conception des systèmes d'onde lumineuse à base de 
soliton est la largeur de l’impulsion T 0 . Chaque impulsion solitonique occupe seulement une 
fraction « temps-bit » et dans la pratique, le lancement de deux impulsions dans deux « temps- 
bit » consécutifs perturbe l'équilibre délicat entre la déformation chromatique et la déformation 
non-linéaire et ceci produit un champ optique combiné qui ne sera pas une solution de l'équation 
non linéaire de Schrodinger (NLS). Ce phénomène est désigné sous le nom de l’interaction des 
solitons adjacents [11]. L’interaction des solitons dans une fibre optique peut être étudiée en 
résolvant numériquement l'équation non linéaire de Schrôdinger (NLS) (III.33) où l’amplitude de 
la paire de soliton à l’entrée de la fibre est décrite sous la forme suivante [11] : 


it(0, r) = sech(r — q 0 ) + p sech(p(r + q 0 )) exp(i0) (III.38) 

T 

Où t = —, p représente l’amplitude relative, 6 la différence de phase initiale et 2q 0 la séparation 
T o 

initiale entre les deux solitons. 
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III.5.2.1. Interaction entre deux solitons voisins d’ordre 1 

Pour étudier l’effet de l’interaction entre deux solitons fondamentaux voisins dans une fibre 
optique sans perte, nous avons utilisé la méthode Split-Sîep-Fourier (SSFM) avec les mêmes 
paramètres de simulations données précédemment. Les résultats de simulations sont présentés 
par la figure (III.6) qui montre l’évolution d’une paire de soliton d’ordre 1 sur une distance z = 
90 L d , avec qo=3.5 et différentes valeurs des paramètres p et 6. 



Figure (III.6) : Evolution d’une paire de soliton d’ordre 1 sur une distance z=90L D . (a) q=3.5, 
p=l, 0=0. (b) q=3.5, p=l, 0=7i/4. (c) q=3.5, p=l, 0=7t/2, (d) q=3.5, p=l.l, 0=0. 

Clairement présentée, l’interaction des solitons dépend fortement de tous les deux paramètres de 
l’amplitude relative et de la phase. D’après la figure (III.6.a), on constate que dans le cas où les 
deux solitons sont en phase (0=0) et leurs amplitudes égales (p=l), ils entrent en collision 
périodiquement le long de la fibre. Cette période est appelée distance de collision donnée par : 
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Lcoi = \ l d ex P(<?o) = h s exp(qr 0 ) (III.39) 

L’effondrement périodique des solitons voisins n’étant pas souhaitable d’un point de vue 
pratique, d’où pour limiter cette interaction, il est nécessaire d’espacer fortement les impulsions 
successives de sorte à avoir L col » L r , où L T représente la distance de transmission. 
Pour q 0 — 8 , nous avons L coi — 300 L s , et cette valeur de q 0 conduit à une réduction du débit 
pour n’importe quel système de transmission optique [11]. 

La figure(III.6.b), montre le cas où les deux solitons présentent une différence de phase (0=7t/4). 
On voit qu’après une première étape d’attraction, les deux solitons s’éloignent l’un de l’autre. 

Dans le cas où (0=7t/2), la Figure (III.6.c) illustre des distorsions de propagation très importants 
des deux solitons et leur écart augmente sans cesse. 

Le dernier cas, Figure (III.6.d), montre l’effet d’une différence d’amplitude (p=l.l). Comme les 
deux solitons sont en phase, ils se propagent en parallèle sans jamais entrer en collision, ni 
s’éloigner l’un de l’autre. Donc de tels petits changements de la puissance maximale (10%) ne 
sont pas nocifs pour le maintenir les solitons. Ainsi, cet arrangement est faisable dans la pratique 
et peut être utile pour augmenter la capacité des systèmes de transmission. 

Egalement la propagation des solitons peut être modifiée par d’autres facteurs, comme l’effet de 
la Gigue de Gordon Haus, les effets supérieurs de non-linéarité et le chirp initial de l’impulsion 
d’entrée [11]. 


III.5.2.2. Solitons d’ordre élevés 

La propagation des solitons d’ordre supérieur (N>1) dans une fibre optique peut être modélisée 
par l’équation non-linéaire de Schrôdinger (NLS). Ces solitons sont caractérisés par une 
puissance plus élevée que la puissance du soliton fondamental (P 0s = N 2 P 0 ) , avec une 
amplitude initiale donnée sous la forme suivante : it(0, r) = Nsech(0, t). Par rapport au soliton 
fondamental, le soliton d’ordre supérieur ne se propage pas identiquement à lui-même quel que 
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soit z, mais présente un comportement périodique , ainsi il retrouvera une forme identique après 
une distance de propagation L s appelée période soliton définie par : L s — -L D [11]. 

Avec la même méthode de simulation numérique (SSFM), et avec les mêmes paramètres de 
simulations que précédemment, nous avons présenté dans la figure (III.7), le comportement 
spatial temporel d’un soliton d’ordre 2 et d’un soliton d’ordre 3 sur deux périodes solitons dans 
une fibre optique monomode et sans perte. 




Figure (III.7) : (a) Evolution d’un soliton d’ordre 2, (b) Evolution d’un soliton d’ordre 3. 


D’après la figure (III. 7, a), nous constatons que le soliton d’ordre 2 et après avoir subi une 
compression, peut retrouver périodiquement sa forme initiale toutes les distances z — L s . Mais 
pour le cas d’un soliton d’ordre 3 et d’après la figure (III. 7, b), nous constatons que le soliton 
dans une première phase subit une compression, puis apparait un éclatement de l’impulsion en 

deux impulsions distinctes à z = —, l’évolution inverse se produit ensuite si bien que les deux 

impulsions se rassemblent. Après une distance de propagation L s l’impulsion redevient identique 
à l’impulsion initiale. Cette évolution périodique peut s’interpréter comme un équilibre 
périodique entre forces attractives (de reformulation) et forces répulsives (d’élargissement). 


7 ^ 
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III.5.2.3. Interaction entre deux solitons d’ordre 2 

Pour étudier l’effet de l’interaction entre deux solitons d’ordre 2 dans une fibre optique sans 
perte, nous avons utilisé la méthode de Split-Step-Fourier (SSFM) avec les mêmes paramètres 
de simulations précédentes. 

L’amplitude de l’impulsion d’entrée est donnée sous la fonne : 

u( 0,t) = Nsech(r — q) + Nsech(r + q) (III.40) 

Les résultats de simulations sont présentés dans la figure (III.8) et la figure (III.9). 



Figure (III. 8): (a)Interaction entre deux solitons d’ordre 2 pour q=3, (b) Profil de l’interaction. 


A partir de la figure (III.8.a), on constate que l'interaction entre les deux solitons est très faible 
quand la séparation initiale est égale à 6. Cette faible interaction induit l'apparition d'un troisième 
soliton de faible intensité par rapport aux autres à la fin de la propagation. Ces évolutions sont 
bien montrées dans la figure (III.8.b) qui présente le profil de cette interaction. 

Par contre si la séparation initiale est diminuée à 3, Figure (III.9.a), dans une première phase les 
deux solitons se compressent avec l’apparition d’un troisième soliton de faible intensité. Puis à 
environ 2 périodes ils s’entraînent en forte collision, puis l’impulsion résultante subit un 
éclatement en trois solitons, où l’intensité du soliton du milieu est plus élevée que celle des 
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autres au cotés. Ces évolutions sont bien montrées dans la figure (III.9.b) qui présente le profil de 
cette interaction. 



Figure (III.9): (a) Interaction entre deux solitons d’ordre 2 pour q=1.5, (b) Profil de l’interaction. 


III.6. CONCLUSION 

La troisième partie de ce travail a été consacrée à l’étude théorique de la propagation non- 
linéaire dans le réseau de Bragg fibré, en particulier la propagation des solitons. Dans un 
premier temps nous avons décrit le concept soliton et les différents types de solitons optiques 
pouvant se propager dans les fibres optiques. La formation des solitons Bragg, structures 
localisés au sein de la bande interdite du réseau a été largement étudiée. Ces solitons étant des 
solutions analytiques de l’équation des modes couplés non-linéaires complètement décrite par 
Aceves et Wabnitz, leur propagation a été modélisée par l’équation de Schrôdinger non-linéaire 
comportant les effets linéaires et non-linéaires. Cette équation a été résolue numériquement par 
la méthode Split-Step-Fourier. Dans ce contexte, nous avons traité numériquement l’impact des 
différents effets pouvant affecter la transmission optique par soliton. L’étude des effets de pertes 
et l’interaction entre solitons voisins ont étés aussi présentés. 




75 








Chapitre : III 


Théorie des solitons. 


III.7. REFERENCES 

[1] N.Zabusky, M.D.Kruskal, « interaction of soliton in collisionless plasma », Phys. Rev. Lett. 
15,240(1965). 

[2] W.Chen and D.L.Mills, “Gap solitons and the nonlinear optical response of superlatices” 
Phys. Rev. Lett. 58 (2), 160 (1987). 

[3] J. S. Russell Report on waves “Report of the fourteenth meeting of the british 
association for advancement of science”, Londres, 311-390, Plates XLVII-LVII (1844). 

[4] Glifford. S. Gardner, John M.Greene, Martin D. Kruskal, and Robert M. Miura, 
“Méthode for solving the korteweg-devries équation”, Phys. Rev. Lett. 19, 1095 
(1967). 

[5] Benjamin & Feir “The disintegration of wave trains on deep water” part 1. Theory, 

Journal of Fluid Mechanics 27, 417-430 (1967). 

[6] Cuenda & Sanchez “Soliton dynamics in a simple model of DNA: Importance of the 
sequence-nonlinear dynamics of spatio-temporal selorganization” Barcelona, UC3M 
(2004). 

[7] Ketterle ‘Nobel lecture’: When atoms behave as waves: Bose-einsteincondensation, 
and the atom laser, Review of Modem Physics 74, 1131-1151 (2002). 

[8] V.E.Zakharov and A.B.Shbat, Zh.Eksp.Teor.Fis. (Soviet) 61,118 (1971); JETP 
(Soviet) 34, 34, 62 (1972). 

[9] A.Hasegawa and F.Tappert “transmission of stationary nonlinear optical puises in dispersive 
dielectric fibers .1. Anomalous dispersion”, Appl. Phys. Lett. 23, 142-144 (1973). 

[10] A.Hasegawa and F.Tappert “transmission of stationary nonlinear optical puises in dispersive 
dielectric fibers .II. Normal dispersion”, Appl. Phys. Lett. 23, 171-172 (1973). 

[11] Agrawal G. P., “Nonlinear Fiber Optics”, 3rd ed. Academie Press, San Diego, (2001). 

[12] L.F.Mollenauer, R.H.Stolen, and J.P.Gordon, “Experimental observation of 
picoseconds puise narrowing and solitons in optical fibers”, Phys. Rev. Lett. 45, 
1095 (1980). 

[13] Y.S.Kivshar, G.I.Stegeman “Spatial optical soliton guiding light for future 
applications” Opt.Photon.News, 13, 5, (2002). 

[14] Wise, Di Trapani “the hunt for light bullet spatiotemporal soliton” Optics and 
Photonics News (Spécial issue soliton) 13, 28-32 (2002). 

[15] Azouzi.Faical .”Etude de la dynamique de la propagation non linéaire des solitons 
optiques modélisés par l’équation de Schrodinguer non linéaire” These de doctorate 
Annaba (2010). 

[16] McLeod et al. « (3+1) -dimensional optical soliton dragging logic », Physical Review A 
52, 3254-3278 (1995). 

[17] C.m.de Sterke and j.E.Sipe. «Envelope function approach for the electrodynamics of 
nonlinear periodic structures ».Phys.Rev.A 38, 5149 (1988). 

[18] J.E.Sipe et H.G.Winful.”Nonlinear Schrodinger soliton in a periodic structure ». Opt. Lett, 
13(2), 132(1988). 





Chapitre : III 


Théorie des solitons. 


[19] A. B. Aceves et S. Wabinitz. « Self-induced transparency solitons in nonlinear refractive 
periodic media ». Phy. Lett. A 141 (1,2), 37 (1989). 

[20] D. N. Christodoulides et R. I. Joseph. « Slow Bragg solitons in nonlinear periodic 
structures ». Phys. Rev. Lett. 62 (15), 1746 (1989). 

[21] J. Feng et F. K. Kenubuhl. « Soliton in a periodic structure with Kerr nonlinearity “. IEEE J. 
Quantum Electron. 29(2), 590 (1993). 

[22] D. Tavemer, N. G. Broderick, D. J. Richardson, M. Ibsen, and R. 1. Laming, “All-aptical 
and gâte based on couped gap-soliton formation in a fiber Bragg grating” Opt. Lett. vol. 23, 
pp. 259-261, (1998). 

[23] B.J. Eggleton, C. M artijn de sterke and R.E.Slusher, J. Opt. Soc. Am. B 16, 587 (1999). 

[24] Agrawal G. P,” Application of nonlinear Fiber Optics”. Academie Press, San Diego, (2001). 

[25] D.J.Kaup and A.C.Newell,”on the coleman correspondence and the solution of the massive 
thirring model” Lettere al Nuuovo Cimento 20,325-331 (1977). 

[26] A.C.Scott, F.Y.Chu, and D.W.Mclaughlin “the soliton: a new concept in applied 
sience”, Proceedings of the IEEE61, 1443-1483(1973). 

[27] Z .B. Wang, H. Y. Yang and Z .Q Li «The numerical analysis of soliton propagation with 
split-step fourier transform method » journal of Physics: Conferenc Sériés 48 878-882 
(2006). 


77 



Chapitre IV : Propagation des soldons Bragg et non-linéarité élevée 


Chapitre : IV 

PROPAGATION DES SOLITONS BRAGG ET NON- 

LINEARITE ELEVEE 


78 









Chapitre IV : Propagation des solitons Bragg et non-linéarité élevée 


IV.l. INTRODUCTION 

L’existence et la formation des solitons dans un réseau de Bragg fibré sont comparables à 
ceux des solitons dans une fibre optique. Ces solitons résultent de deux effets combinés qui 
sont l’effet Kerr (auto-modulation de phase) et la dispersion chromatique induite par le 
réseau. 

Le soliton Bragg peut se propager sans déformation dans le réseau avec des fréquences 
proches ou appartenant à la bande interdite mais pas en dehors de celle-ci. Dans la littérature 
scientifique, les solitons Bragg sont appelés (Gap soliton) quand leur fréquences sont situées 
dans la bande interdite. Plusieurs études expérimentales ont démontré l’existence des solitons 
Bragg avec une caractéristique principale, celle de pouvoir se propager à des vitesses 
largement inférieures à celle de la lumière, des vitesses inférieures de près de 50% ainsi ont 
été démontrées [1-5]. 

Théoriquement, la propagation des impulsions optiques dans un réseau de Bragg est 
modélisée par les équations des modes couplées non-linéaires (NLCME) [6] et les solitons 
Bragg sont les solutions analytiques de ces équations. Christodoulides et Joseph sont les 
premiers à avoir développé une méthode basée sur le modèle de Thirring massif pour 
rechercher des solutions situées exactement au centre de la bande interdite [7]. Peu après, 
Aceves et Wabnitz [1] ont généralisé cette méthode afin de déterminer des solutions pour 
toute la bande. 

Récemment, plusieurs études théoriques et numériques ont été réalisées pour clarifier et 
classifier les conditions de génération des solitons Bragg, et la dynamique non-linéaire de 
leur propagation [8-10]. 

Dans ce chapitre nous avons étudié analytiquement l’existence des ondes solitaires dans un 
réseau de Bragg fibré en régime femto-secondes, où la propagation des impulsions est 
modélisée par un système d’équations des modes couplés, présentant des termes d’ordre 
élevé, des non linéarités de troisième, cinquième et septième ordre. Les applications 
physiques de ce modèle se situent dans un domaine de pointe qui est celui de la génération 
des solitons Bragg dans un milieu fortement non linéaire comme le verre de chalcogénure. Ce 
type de verre présente non seulement des non-linéarités de troisième et cinquième ordre mais 
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aussi de septième ordre [11,12]. La résolution de ce système d’équations non intégrables a 
nécessité l’application de la méthode perturbative multi-échelle[13,14] qui a permis la 
réduction du système d’équations non linéaire couplées (NLCME) en une seule équation dite 
équation de Schrôdinger non linéaire perturbée (PNLS) qui décrit la propagation non linéaire 
des impulsions aux bords de la zone du gap de la bande photonique. Vu que cette dernière 
équation présentait des termes non-Kerr et des termes de dispersion de troisième et quatrième 
ordre, sa résolution a nécessité l’utilisation de la fonction d’essai de Li et al. [15]. 

Les solutions exactes obtenues sont des solitons noirs et brillants combinés. Aussi, nous avons 
également présenté les conditions d’existences de ces solutions ainsi que leurs importants 
paramètres physiques tels que la puissance et la largeur de l’impulsion. 


IV.2. EQUATIONS DES MODES COUPLES 


La variation de l’indice de réfraction dans un réseau de Bragg en verre de chalcogénure est 
donnée par cette relation [12] : 


n ( z) — n +Hj( z )cos[ — + z + 0( z )] + h 2 \e(z)\ 2 +h 4 \e(z)\ 4 +n 6 \E{zf 


(IV.l) 


m, n 4 et n6 représentent les coefficients de la non-linéarité avec m et n6 > 0 et n 4 <0. h 
représente la modulation moyenne de l’indice de réfraction, n x (z) l’amplitude de modulation 
de l’indice le long de la fibre, A la période de la modulation et 0(z) sa phase. L’évolution du 
champ électrique E(z) dans le réseau de Bragg fibré prend la forme suivante [16]: 


E(z, t ) = [A + (z, t)exp(ik B z) + A_(z, t)exp(- /fc B z)]exp(- ico 0 t) 

(IV. 2) 

Où A+ et A- sont les amplitudes de l’onde progressive et de l’onde régressive respectivement. 
Jc B est le vecteur d’onde de Bragg et coo la fréquence centrale. 

Pour établir le système d’équations des modes couplés non-linéaires d’ordre élevé décrivant la 
propagation des impulsions femto-secondes dans le réseau de Bragg, on introduit les 
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équations (IV. 1) et (IV.2) dans l’équation de Helmotz (III. 10), obtenue à partir des équations 
de Maxwell. Après quelques simplifications et quelques transformations algébriques on 
obtient le résultat suivant [12]: 


f ( 1 8A + (z,t) + 5A + (z,tŸ 


V v s dt 


dz 


■/c4_(z,0 + r i (A + | 2 +2|A_| 2 )A + 


y 


r 2 [|A + | 4 + 6 |a + | 2 |a_| 2 +3|a_| 4 )a + +r 3 |A + | 6 +12 |a + | 4 |a_| 2 +18|a + | 2 |a_| 4 + 4|aJ 6 1a^ =0 


5 k 


^ 1 dA_(z.,t ) dA_(z,t)' 


V v * dt 


dz 


+ k4 + (z, t) + T, j A_ | 2 + 2|.A + 1 2 )a_ 


+ 


y 


(IV.3) 


rdl A_ | 4 + 6| A_ | 2 1A + 1 2 + 3| A + 1 4 )a_ + r 3 1 A_ | 6 +12| A_ | 4 1A + 1 2 +18| A_ | 2 1A + 1 4 + 4| A_ | s )a_ = 0 


Le système d’équations (IV.3) modélise la propagation des impulsions optiques ultras courtes, 
dans un réseau de Bragg fibré ou t et z représentent le temps et la distance normalisés 
respectivement, k est le coefficient de couplage linéaire et v g la vitesse de groupe. Les termes 
dépendant de T] représentent la non-linéarité cubique (la modulation des phases croisées et 
l’auto modulation de phases), les termes dépendant du paramètre T2 décrivent la non-linéarité 
de cinquième ordre et les termes dépendant du paramètre T 3 représentent la non-linéarité de 
septième ordre. 

Dans le cas linéaire c’est à dire (F]=0, r2=0 et T3=0), le système d’équations (IV.3) peut être 
résolu facilement et la solution prend la forme suivante [14,17]: 



U exp (i(kz~ cot))+ c.c, U 


U + 

\U-J 


(IV .4) 


Où le vecteur d’onde k et le déphasage de fréquence co satisferont la relation de dispersion 
linéaire suivante: œ 2 — k 2 + k 2 . Dans le cas particulier où k = 0, œ — 1 k et la solution 
(IV.4) sera comme suit: 


X 

U- 


y 


= c 


f i 3 

exp(-?'/rf)+cc 

V V 


(IV.5) 
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La solution (IV.5) satisfait à: 



Où L est un opérateur donné par: 



k id, 


(IV .6) 


(IV.7) 


En effet la résolution des équations des modes couplées non linéaires (IV.3) s’avère difficile 
du fait que ces équations ne sont pas intégrables au sens mathématique. Cependant dans la 
limite des faibles intensités, ces équations peuvent se réduire en une seule équation, celle de 
Schrodinger non-linéaire perturbée (PNLS), par l’application de la méthode perturbative 
multi-échelle [14, 17-21]. Cette équation étant intégrable et ses solutions sont de véritables 
solitons [17]. 


IV.3. EQUATION DE SCHRÔDINGER NON-LINEAIRE PERTURBEE (PNLS) 

Les solutions recherchées dans le cas non-linéaire du système d’équations des modes couplés 
(IV.3) en appliquant la méthode perturbative multi-échelle sont de la forme [8, 14,17]: 


(A ) 


m 

+ 

| = e 1/2 A(r 1 ,r 2 ,Z)| 




K 


exp(- iKt) + sU x + 


S 3/2 U , 


+ £ 2 U , 


+ ... 


(IV .7) 


Où xi = et, i 2 =e 2 t, Z= £ 1//2 z. 


La résolution du système (A+, A.) se fera suivant les ordres croissants de s et le 
développement en fonction du terme O(s) donne: 




vly 


LU, =~id z A | 

La solution de ce problème linéaire est donnée par: 


(IV.8) 


U , =- d.A 

2 k 1 


f L 
\b 


(IV.9) 
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Nous continuons notre démonstration pour les corrections d’ordre élevées de (A+, A.) et la 
solution (IV.9) sera utilisée pour le calcul des termes non-linéaires: 

> Le développement des termes représentant la non-linéarité cubique donne: 


[|A 2 + 2A | 2 k = 

e'/'-A-i— — 

2 

+ 2 

V 

2 k dz 



— e 


l/2 A-i —— 

2 K ôz 


2 \ 


r 1/2 A-i —— 

2 K dz 


-ÏKt . 

e + c.c. 


{a | 2 + 2|aJ 2 )a = 

. e <r- A _ t ^ôA 

2 

+ 2 


2 ^ 

('- s V2 A _(2L®f| 


M -1 l +1 r - 

V 

2k dz 


2k dZ 

/ 

(_ 2k dZ z 


En simplifiant les termes des modules aux carrés, on obtient : 

(|A,| 2 +2|A_| 2 )a, =[(3^ 2,2 |A| 2 A-/|-(2|A| 2 a z A + A 2 8 2 A , )+ j^{lAd z Aa z A‘ + A'(d z Af) 
-i^(d z A)%A')y i " + c.c 

(|A_| 2 + 2|A^|’)a_ =-[(3i 3,2 |A| ! A+i|-(2|A| 2 S z A + A 2 0 z A , )+|- T (2A0 z Aa z A‘ + A"(0 z A) 2 ) 

+i£j-(a z A) 2 (a z A-)] e -'“+^ 


> Le développement des termes présentant la non-linéarité de cinquième ordre donne: 

(a$ + 6A j ./ ! A ' + 3A 4 )/t. = [l0f 5/2 |A| 4 A-/ — (6\A\ 4 d z A + 4A 2 \A\ 2 d z A*) 

2/c 

+ ^(2A 3 (ô z A*) 2 +6A*|A| 2 (a z A) 2 + 12A|A| 2 (ô z AXô z A*)) 

- i ^ (ôA 2 (d z A)(ô z A* ) 2 + 12 |a| 2 (d z Àf (ô z A* )+ 2A* 2 (d z A)(ô z A) 3 ) 

+ ^(6A(Ô Z A) 2 (<3 Z A* ) 2 + 4A* (d z Àf{d z A* )) 

~ i -^ T (lO(ô z A) 3 (ô z A*) 2 ) +c.c. 
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|a_| 4 + 6|a_| 2 |a + | 2 +3|a + | 4 )a_ = 

[-l0s 5,2 \A\ 4 A-i^-(6\A\ 4 d z A + 4A 2 \A\ 2 d z A , ‘) 

-^(2 A 3 (d z A*) 2 +6A*|A| 2 (5 z A) 2 +12A|A| 2 (5 z A)(ô z A*)) 

-/^(ôA 2 (5 z A)(ô z A*) 2 + 12|A| 2 (5 z A) 2 (5 z A*) + 2A* 2 (5 z A)(ô z A) 3 ) 

oK 

-^(6A(5 z A) 2 (ô z A*) 2 +4A*(ô z A) 3 (5 z A*)) 

-'■3 ^(lO(d z Aj(d z A*) 2 )y“ +C . C . 


Le développement des termes présentant la non-linéarité de septième ordre donne: 

|a + | 6 +12|a + | 4 |a_| 2 +18|a + | 2 |a_| 4 +4|a_| 6 )a + = 

[35s 7/2 |A| 6 A - i y- ( 3 A 2 1 A| 4 5 Z A* + 4| A| 6 5 Z a) 

+ ^(a 3 |A| 2 (ô z A*) 2 +2A*|A| 4 (5 z A) 2 +4A|A| 4 (5 z AXd z A*)) 

- i (l2A 2 |A| 2 (5 Z AXô z A* f +18|A| 4 (5 z A) 2 (d z A* )+ 4A* 2 1A| 2 (d z A) 3 + A 4 (d z A* )'’) 


(l8A|A| 2 (5 z A) 2 (5 z A* j 2 + 12A*|A| 2 (5 z A) 3 (d z A* )+ 4A 3 5 z a(<3 z A* J + A* 3 (d z A) 4 ) 


1 6at 

,15g 6 

32k 5 


(4|A| 2 (5 Z A) 3 (5 Z A*) 2 + 2A 2 (S Z A) 2 (<3 Z A*) 3 + A' 2 (ô z A*)(ô z A) 4 ) 
^(4A(5 z A) 3 (5 z A*) 5 +3A*(S z A) 4 (d z A*) 2 ) -‘^^((Ml^A*) 3 ) K 


128/f ' 


+ C.C. 


|a_ | 6 +12| A_ | 4 1A + 1 2 +18| A_ | 2 1A + 1 4 + 4| A + 1 6 )a_ = 

[- 35f 7/2 |A| 6 A - i — (3A 2 |A| 4 ô z A* + 4|A| 6 ô z A 
2k 


I5e 9/2 

4k 2 
,3g 5 
1 8k 3 


(a 3 |A| 2 (ô z A*) 2 +2A*|A| 4 (5 z A) 2 +4A|A| 4 (a z A)(ô z A*)) 
(l2A 2 |A| 2 (ô z A)(ô z A*) 2 +18|A| 4 (ô z A) 2 (ô z A*)+4A* 2 |A| 2 (3 z A) 3 + A 4 (d z A*) 3 ) 
^(l8A|A| 2 (ô z A) 2 (ô z A*) 2 + 12A*|A| 2 (5 z A) 3 (5 z A‘)+4A 3 ô z A(ô z A*) ? + A* 3 (5 z A) 4 ) 
■(4|A| 2 (ô z a) 3 (ô z A*) 2 +2A 2 (ô z A) 2 (a z A*) 3 + A* 2 (ô z A*](ô z A) 4 ) 

( 4 A(a z a) 3 (ô z .a*) 3 + 3A*(a z a) 4 (d z A*)") -i ((d z A) 4 (<9 z A*) 3 ) ]e 


16* 

,15g f 


32k 

5g‘ y2 
" 64k 6 


+ C.C. 
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La balance des termes de 0(£ 3 ^ 2 ) donne : 


LU 2 = 


( 1 , A 

-id. A -ô 2 2 A-3r,U|"A 

^ ri 2 k z 11 


f i 3 


A-V 


i —ixt i 

e + cc. 


(IV. 10) 


Il est à noter que pour une faible variation des termes autour de U 2 et avec (_* )e lKt dans 
l’espace non nul de L, nous obtenons la condition physique suivante: 


id T A +—d\A + 3Y ] \A\ 2 A = 0 
Tl 2 k z 11 1 


(IV. 11) 


L’équation (IV. 11) représente l’équation de Schrôdinger non-linéaire (NLS), qui décrit la 
propagation des impulsions en dehors de la bande interdite du réseau. Cette équation est 
complètement intégrable par la méthode de la diffusion inverse. 

Dans le cas où les effets d’ordre élevé seront considérés, les termes perturbatifs présents dans 
l’équation (IV.3) doivent être ajoutés à l’équation de Schrôdinger (NLS). 


En continuant cette étude sur l’impact des effets non-linéaires d’ordre élevé ainsi que la 
balance des termes 0(s ) et en utilisant l’équation (IV. 11), on trouve : 

f / \ i V | N 


lu 3 = 


V 2/c 

La solution de cette équation est : 

1 


i ^ (21 Al 2 d. A + A 2 d . A* )+ Ar A 

' , - v 1 1 2 ‘ 7 4 k 2 z 


vly 


. ~lKt , 

e + cc. 


U : , = -i 


4k 2 


2^(1 A\ 2 A) + —d\A 

2 AT 




ai y 


\e ,Kt +cc 


(IV. 12) 


(IV. 13) 


Par conséquence l’équation non-linéaire de Schrôdinger perturbé (PNLS) résultante s’écrit 
comme suit : 

id r A + — d 2 2 A + 3r,\A\ 2 A + i-^—d\A + i^L-d(\A\ 2 A)=0 (IV.14) 

2 k z 11 1 8/c 3 z 2 k 2 ^ J 

On note que ce modèle d’équation a été déjà traité par K. Mezghiche et al. [8]. Ils ont adopté 

la méthode de l’amplitude complexe des solitons combinés présentée par Li et al. [15] et les 

solutions trouvées sont des ondes solitaires analytiques de type soliton brillant et soliton 

sombre. 
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Pour mieux décrire la génération des solitons Bragg auprès de la bande photonique (PBG) 
dans un milieu périodique a non linéarité Kerr, nous continuons le développement des 
termes 0(£ 5 / 2 ) . 

L’équation résultante se présente comme suit: 

id T A + — d\A + 3r.\A\ 2 A + i-^—d\A + i-^-Ô Ja| 2 a)+—*— 5 4 4 A + ^-Ma| 4 A = 0 
r ' 2k z 11 8a- 3 z 2 K 1 A 1 ’ 16a- 4 z k 

(IV. 15) 


Cette équation a été résolue analytiquement par Senthilnathan et al. [17]. la méthode de 
couplage amplitude-phase a été adoptée et la formation d’un soliton brillant a été démontrée. 


T 

En continuant l’ajout des termes perturbatifs 0(s ), nous arrivons à l’équation résultante 
suivante: 


id T A + — — ô 3 ,A + 3r. |a| A + i 
T ' 2 k z 11 1 


-ô 3 A + i 
8a- 3 z 


—-rd \tA\ 
2k~ z 


A + 


-î—r d 4 4 A + 

16a* 4 z 


5F, 


K 


AA 


-i 


2k 1 


5.1 A A =0 


(IV. 16) 


Et finalement, toujours en adoptant la même méthode concernant les termes 0(£ 7//2 ) , on 
obtient l’équation non-linéaire de Schrodinger perturbé (PNLS) a non-linéarité d’ordres 
élevés suivante: 

iô_A+a x d\ A+a 2 |A|" A+ia 3 ô 3 t , A + afî^ A+ia 5 d t |a|" a|+û: 6 |A| 4 A — cr 7 5, (| A| 4 A j+ 
a s | A| 6 A - ia 9 d t | A| 6 a|= 0 

(IV. 17) 


L’équation (IV. 17) décrit la propagation non-linéaire des impulsions femto-secondes voisines 
mais en dehors de la bande interdite du réseau de Bragg fibré où ti = z et Z= t. Les 
paramètres physiques a, (avec i=l,..,9 ) représentent des paramètres réels donnés par les 
relations suivantes : 


ai “2 K ’“ 2-3ri ' a3 “8K3' a4_ 16K« 

r 2 35IN 5 IL 

a 7 - —-, a 8 =-, do = — 

' 2k 2 a 2 k y 4 k 2 



«6 


sr 2 


t 

K 


(V.18) 
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L’équation (IV. 17) inclut les effets non-linéaires et les effets de dispersion d’ordres 
supérieurs. Les termes proportionnels a a lt a 3 a 4 sont la dispersion (GVD), le troisième ordre 
de dispersion (TOD) et le quatrième ordre de dispersion (FOD) respectivement. Les termes 
proportionnels a a 2 ,a 5 sont l’auto-modulation de phase et l’auto- raidissement de phase 
respectivement. Les termes proportionnels a a 6 , a 7 , représentent les effets non-linéaires de 
cinquième ordre. Les termes proportionnels a a 8 ,a 9 , représentent les effets non-linéaires de 
septième ordre. 

Ce nouveau modèle d’équations (IV. 17), que nous avons atteint présente notre base de travail 
de recherche. Cependant ce modèle d’équation n’est pas complètement intégrable à cause de 
l’ordre supérieur de la non-linéarité et de la dispersion en même temps, et ne peut pas être 
résolue par la méthode de la transformation inverse. 

Pour résoudre ce modèle d’équation, on a adopté la méthode de l’amplitude complexe des 
solitons combinés présentée par Li et al. [15]. Cette méthode permet de dériver la forme de la 
solution combinée qui est la superposition entre deux amplitudes de soliton noir et brillant, de 
tel sorte que la forme W-Shape [8,15] de l’onde solitaire apparait avec des conditions sur les 
paramètres physiques de l’équation non-linéaire. 


IV.4. RESOLUTION DE L’EQUATION NON-LINEAIRE PERTURBEE (PNLS) 

Afin d’établir les solutions analytiques des ondes solitaires de l’équation de Schrôdinger non- 
linéaire perturbée (IV. 17) que nous avons obtenu précédemment, la méthode de l’amplitude 
complexe des solitons combinés présentée par Li et al. [15] a été utilisée. Cette méthode 
présente la forme combinée de la solution qui est la superposition entre deux amplitudes de 
soliton noirs et brillants, avec une valeur asymptotique ne s’approche pas au zéro lorsque le 
temps s’approche à l’infini. 

La solution de l’amplitude complexe s’écrit sous la forme [15] : 

A(z,t) = q(z,t)exp(i0 L (z,t)) (IV. 19) 
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Où q(z,t) présente la fonction de l’enveloppe complexe et </> L (z,t) = kz-cot le shift linéaire de la 
phase. En insérant cette solution dans l’équation (IV. 17) et en enlevant le terme exponentiel, 
l’équation résultante prend la forme suivante : 


id ,q + ci x q - ia 2 ô t q + a 3 d\q + ia 4 d 3 ,q + a 5 \qf q + ia 6 \q\~d r q + ia 1 q 2 d l q* + « g |g| 4 g 
-ia 9 \q\ 4 ô t q - ia l0 q 2 \q\~ d t q* +a n \q\ 6 q-ia 12 \q\ 6 d t q-ia l3 \q\ 4 q 2 d t q* + a l4 d 4 4 q = 0 


(IV. 20) 


Avec : 

a { = -{k + a i m 2 + a 3 co 3 + a 4 m 4 \a 2 = (2 a l co + 3 a 3 af -4a 4 m 3 \a 3 = (a { + 3 a 3 m- 6 a 4 co 2 ), 
a 4 = a 3 - 4a 4 co,a 5 = ( a 2 + a 5 co),a 6 = 2 a 5 ,a 2 =a 5 ,a 8 = (a b - a 1 co),a g = 3 a 7 ,a l0 = 2 a 7 , 

1 = ( OCÿ OCgÛ)), Cl^ 2 — 4 OCq , CZj j = 3 CCg ,cty 4 = ct 4 

(IV.21) 

Nous recherchons mantenant les solutions des ondes solitaires, en utilisant la fonction d’éssai 
(Ansatz) présentée par Li et al suivante [15] : 

q(z,t) = i/3 + A tanh(P 0 (t-/3 x z)) + ipsec h[P Q (t - /3 x z)\ (IV.22) 

Où les paramètres Po et fii sont la largeur de l’impulsion et le décalage de l’inverse de la 
vitesse de groupe respectivement. Ces paramètres prennent des valeurs réelles, mais /?, p et A 
peuvent être des réels ou des nombres complexes. Il est important de mentionné que l’ansatz 
ci-dessus a également été appliqué avec succès pour certain modèles d’équations (NLS) 
contenant des termes d’ordre élevé de la non linéarité et de la dispersion [22-24]. 

Cependant l’amplitude de l’impulsion s’écrit comme suit : 

\q{z,t\ ={â 2 + j3 2 )+2j3psech[P Q (t-J3 { z)] + (p 2 -A 2 )sec/i 2 [P 0 (f-/?,z)] (IV.23) 

Et la fonction du déphasage non-linéaire est donnée par : 


f 

<fi NL (z,t) = arctan 

V 


P + psech[P 0 (t- Pp)\ 

2 tan h [/),(/ -//,-)] , 


(IV. 24) 


Nous commençons notre méthode de résolution par la substitution de l’équation (IV.22) dans 
l’équation (IV.20). En séparant les termes de tanh des termes de sech, des coefficients aux 


88 




Chapitre IV : Propagation des solitons Bragg et non-linéarité élevée 


termes indépendants contenants des combinaisons indépendantes de ces fonctions 
hyperboliques apparaissent et sont égales à zéro. 

Après séparation des parties réelles et imaginaires, nous obtenons les 17 équations 

paramétriques suivantes : 


+ afA + /? 2 ) + a 8 (/f + p 2 ) 2 + a n 


[X-+P 1 -) 



M 


p{a t + a 2 P 2 + ci 5 [a 2 +3/3 2 )-2a 7 /3AP 0 + a 8 (/l 2 + /3 2 )(a 2 +5/3 2 )+2a l0 (A 2 + /3 2 ]A/3P 0 +a n (A 2 + /3 2 J(a 2 +7 /? 2 )+ 

2a n (A 2 +/3 2 ) 2 A/3P 0 +a H P 0 4 }=0 (2a) 


PA - aA 0 + 4 a A P 0 ' + a 6 {A + P 2 )p o + 4 i 2 - 2p 2 - p 2 )p o - a 9 (A + p 2 )P 0 - ajA - 2p 2 - p 2 \X~ + p 2 )-4p 2 p 2 ]p a 
-a n (A 2 + p 2 7P 0 -a n [(A 2 -2p 2 -p 2 )(A 2 + p 2 )-Zp 2 p 2 ]p o {a 2 + p 2 )} + p{a 5 [3p 2 -A^+a^A 2 + p 2 ){3p 2 - A 2 ) + 4p 2 p 2 \ 


+ 3a n {A 2 + p 2 lA 2 + p 2 ){3p 2 -A 2 )+4p 2 P 2 ] } = 0 


(3a) 


p{-2a 3 P 2 + a 5 (p 2 -X-)+2a 6 ApP 0 +a s \2{p 2 - A 2 \A + 3/7 2 )+4/?V]-4a 9 (/l 2 + p 2 )ApP ü 
- 2a w (2A 2 -3 p 2 - P 2 )APP 0 + a n [3{p 2 - A 2 \A 2 + p 2 \A + 5 p 2 )+4(3 A + 5 p 2 )p 2 p 2 ] 

-6a l2 (A 2 + P 2 ) 2 APP 0 -4a 13 [(/t 2 + p 2 j2A 2 -3p 2 -p 2 )-2p 2 p 2 ]ipP 0 -20a u p 4 0 } = 0 (4a) 


A{-6a 4 P 0 i +{a 6+ a 1 ip 2 -A 2 )P 0 -2a 9 [{p 2 - A\A + p 2 )+2p 2 p 2 ]p o + 2a lo {p 2 -A 2 Jp o 
- 3a l2 {A 2 +p 2 l{p 2 -A 2 lA 2 +p 2 )+4p 2 p 2 ]p 0 -a43A 2 -4p 2 -p 2 l(A 2 +p 2 lp 2 -A 2 )+ 4p 2 p 2 ]p o } 

+ p{a i {p 2 -A 2 \5p 1 -A 2 )+a l fp 2 -A\ (a 2 + p 2 \np 2 -3A 2 )+4p 2 (A + 4p 2 ) }+8 p 2 p 4 ] } = 0 (5a) 


p{ a s (p 2 -A 2 ) 2 -4a 9 (p 2 -A 2 )ApP 0 -2a w (p 2 -A 2 )ApP 0 + 3a n (p 2 -Afp 2 - A\A+3p 2 )+4p 2 p 2 


-4a u [3(p 2 -A 2 \A 2 + p 2 ]+2p 2 p 2 ]pAP 0 + lOa l3 (p 2 - A f AfiP {) + 24a, 4 P 4 } = 0 (6a) 

A{ —a 9 (p 2 -A 2 ) 2 P 0 -a l0 {p 2 -A 2 ) 2 P 0 -3 a l2 {p 2 -Afp 2 -A\A +/? 2 )+4/?V]p h 
- a 13 {p 2 -A 2 fp 2 -A 2 \3A 2 -2p 2 +P 2 f4p 2 p 2 ]p () }+Pa n {p 2 -A 2 J{lp 2 -A 2 )= 0 (7a) 

p{ a n [p 2 -/l 2 ) 3 -6 a n (p 2 - A 2 ) 2 APP 0 -4 a 13 (p 2 - A 2 f A[iP {) } = 0 (8a) 

~ a l2 {p 2 — A ) 3 P 0 — a l3 (p 2 - A ) 3 P 0 } = 0 (9a) 
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A {«! + a,( â 2 + 0 2 )+ a s (A 2 + 0 2 ) 2 + a n (A 2 + 0 2 ))= 0 (lOa) 


p{-PA-a 2 P 0 +a 4 P 3 +2 a 5 A0+a 6 (A 2 +0 2 )p o -a 7 (A 2 -0 2 )p n + 4a t (A 2 + p 2 ]Ap~a g (A 2 + p 2 f P n +a w (A 4 - 0 4 )p n + 
6a n (A 2 + P 2 )Ap-a n (X- + 0 2 }p o +a n (s l 2 + 0 2 \a 4 - 0 4 )p ü }= 0 (lia) 


A{-2a,P 2 +a,{p 2 -X-)+2a,\X- + P 2 \p 2 -0)+ip 2 p^+ïa^X- + P 2 \X; + P 2 \p 2 -X:)+Ap 2 p 2 )-UJ 
+ p{2a 6 p 2 P 0 - 2a-, {a 2 - p 2 )p o - 4a, (a 2 +p 2 )p 2 P 0 + 2 ajx 2 + 0 2 \a 2 - p 2 )+ p 2 {a 2 -0 2 )]p 0 - 6 ü u (a 2 + p 2 J p 2 P 0 


2a, 


Ia 2 +P 2 ) 2 (a 2 -p 2 )+2p 2 (a 4 -P 4 )Ip o } = 0 


(12 a) 


p{-6a 4 P 0 3 + (a 6 + a 7 (p 2 - A 2 )p o + 4 a s (p 2 - A 2 ]a0 - 2a, [{a 2 + 0 2 \p 2 - A 2 ) +2 p 2 0 2 ]p o 
+ 6a l0 (/l 2 -p 2 )p 2 P 0 + 4a^{p 2 -A 2 \a 2 + 0 2 \r 20 1 p 2 \t0 - 3a 12 (/L 2 + 0 2 \a 2 + 0 2 \p 2 - A 2 ) + 4 p 2 p 2 ]p„ 

+ a l3 [l0^ 2 (A 2 - P 2 \a 2 +P 2 )+4P 2 P 2 (a 2 -P 2 )+{p 2 -A 2 \A 4 - 0 4 )\i ü }=0 (l3a) 


A{a s (p 2 -A 2 ) 2 +3a n (p 2 -A 2 l{p 2 - A:\a 2 + 0 2 )+40 2 p\24a l4 P Q 4 }+p{-4a g (p 2 ~A 2 )p 2 P 0 
+ 2 a 10 (A 2 -p 2 \2p 2 -A 2 )p 0 ~4a 12 [3(p 2 -A 2 \a 2 + 0 2 )+ 20 2 p 2 ]p 2 P o 

+ 2a n (X- - p 2 %A: + 0 2 )(3p 2 -2A:)+ p 2 {p0 2 - } = 0 (14a) 

p{-a 9 (/ o 2 - /l 2 ) 2 P 0 - a 10 (/l 2 - p 2 ) 2 P 0 + 6 a u (p 2 - /l 2 f Ap - 3 a 12 (p 2 - A 2 f p 2 - /l 2 \a 2 + 0 2 )+40 2 p 2 ]/), 

+ a 13 (/l 2 - p 2 J [p 2 - A: X^ 2 + P 2 2p 2 (lp 2 - 5A 2 )]p 0 } = O (l5a) 


A{ a n (p 2 -A 2 J }+p{-6a l2 (p 2 -A 2 Jp 2 P 0 -2a 13 (A 2 -p 2 ) 2 (3p 2 -A 2 )P 0 } = 0 


(l6a) 


p{ -a l2 {p 2 -A 2 JP 0 +a l3 (A 2 -p 2 JP 0 } = 0 (17a) 

A l’aide de ces équations paramétriques, nous allons discuter les différents cas possibles 
nous permettant de distinguer les types de solitons, solution de l’équation (PNLS) d’ordre 
élevé introduisant les termes non-linéaires de troisième, cinquième et septième ordre. 
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IV.5. SOLUTIONS POSSIBLES DE SOLITONS BRAGG 

Pour obtenir les solutions des ondes solitaires de l’équation non-linéaire de Schrôdinger 
perturbé (IV. 17), il faut imposer quelques restrictions au paramètres (p, X), de tel sorte que 
les équations (la)-(17a) deviennent compatibles. Pour cela, nous avons le cas des équations 
(7a)-(9a) et (14a-17a) qui sont compatibles pour: p=±L , avec et (a 12 + a 13 ) ^ 0. 
Si nous choisissons le cas p=k, la solution obtenue prend la forme suivante : 

q(z,t ) = i(j + ltanh[ P 0 (t - (pz)] + ilsech[P 0 ( t - p,z )] (IV.25) 

L’un des paramètres physique de cette solution, la largeur de l’impulsion P 0 peut être 
déterminée à partir de l’équation (6a) ou (14a) d’où on trouve : 


P 0 = 


f \P 

a. 


\P a u j 


pl 


Avec l’utilisation des équations (IV.18) et (IV.21), l’équation (IV.26) devient : 

3 

P'y l 


Po = 


r sor 3 K 2 ^ yi 


v ^ J 


(IV.26) 


(IV. 27) 


Le vecteur d’onde k peut être déterminé à partir des équations (la) ou (10a) : 

k — — (cq to + (X> O) + oc 4 co (cx- ; + CX 3 O) )(/r + p — tXjto )(/. + p j + (tx 3 — Otg O) )(/. + fi j 


Avec l’utilisation de l’équation (IV. 18), le vecteur d’onde k s’écrit comme suit : 


m co co 

v 2v 8 V ïëZ, 


1 2 k 1 p ’ 


5L r 2 


co 


y K 2 K j 


(a 2 + p 1 ) 2 


35L 5L 


(IV.28) 


2 k 4 r 


-co 




(IV.29) 


L’équation (1 la) va pennettre l’obtention du paramètre ‘/V: 


/?, = -a 2 + a 4 Pg + 2a 5 P 0 l ip + [à 2 + p 2 \a 6 + 4a g P 0 ' ip + a l3 (l 4 - p A ))+ ( À 2 + p 2 ) (ôu, X P 0 l ip -a 9 ) 
-a\l 2 -p 2 )+a w (X -p A )-a n (^ + P'J 


(IV.30) 
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Et avec l’aide des équations (IV.18) et (IV.21), on trouve : 


A = 


( co 3 3 co 2 co ^ 


v 4 k 4 8 / r 3 k j 

r, ( 2 or 2 2 F 


+ ( 6r '-/ ,: )+E(i‘ -p' )+ 


K \ K 


K 


CO 


pït>i ++/>’)+((^ - 1 - fjJF+yl 


5r, 


a: - 


(i 2 +/? 2 ) 3 


(IV.31) 


Finalement, à partir de l’équation (5a) on obtient ‘co’ : 


a, a 7 
co =—— + —— 
3 a 4 2 ar 9 


4a 9 

V 3°'4 y 


1 

3 3 

+ - 
2 


f 4 a 9 A 
V 3 <2; 4 y 


(4 2 + /E)-^ 

v ’ 3a, 


(IV.32) 


En substituant les coefficients a, (avec i=l,..,9) de l’équation (IV. 18) dans l’équation 
(IV.32) on trouve: 


co = 


8or 3 s- 


2 > 


5-* &**/>■) 


v 5r s 2 


-4 K 


(IV.33) 


Et l’intensité de l’onde solitaire obtenue est donnée par : 
\q(z, t )| = (l 2 + p 2 ) + 2fiX sec h[p 0 (t - fi ; z)] 


(IV.34) 


On constate que la solution obtenue (IV.34) admet deux types de solitons : 

si P A < 0 ,1a solution prend la forme d’un soliton noir, par contre si /?A > 0 la solution prend 
la forme d’un soliton brillant. 


IV.6. CONCLUSION 

Dans ce chapitre, nous avons étudié analytiquement l’existence des ondes solitaires dans un 
réseau de Bragg fibré en régime femto-secondes, où la propagation des impulsions est décrite 
par un système d’équations des modes couplés incluant des termes d’ordre élevé de non 
linéarités de troisième, cinquième et septième ordre. Afin de résoudre ce système d’équations 
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non intégrables, nous avons utilisé en premier la méthode perturbative multi-échelle qui a 
permis la réduction de ce système d’équations non linéaire couplées (NLCME) en une 
équation dite l’équation de Schrôdinger non linéaire perturbée (PNLS). Cette dernière décrit 
la propagation non linéaire des impulsions aux bords de la zone du gap de la bande 
photonique. La résolution de cette équation compliquée contenant des termes non-Kerr et des 
termes de dispersion de troisième et quatrième ordre a nécessité l’application de la fonction de 
l’amplitude complexe combinée. Les solutions obtenues sont de type solitons noirs et 
brillants. En outre, nous avons également présenté les conditions d’existence de ces solutions 
ainsi que leurs importants paramètres physiques tels que la puissance et la largeur de 
l’impulsion. 
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Conclusion générale. 


Les travaux de recherche menés au cours de cette thèse ont majoritairement porté sur une 
étude théorique et numérique de la propagation non-linéaire des impulsions solitoniques dans 
la fibre optique et le réseau de Bragg fibré et ceci dans le cadre d’applications en 
télécommunications optiques à haut débit. 

Notre étude nous a conduites en particulier au traitement de la dynamique non-linéaire des 
solitons Bragg dans un réseau de Bragg fibré en régime femto-seconde, où la propagation des 
impulsions pouvait être modélisée par un système d’équations des modes couplés présentant 
des termes de non linéarité d’ordre élevé. Cependant, l’introduction de ces termes change la 
nature du système qui devient non intégrable. L’obtention des solutions correspondantes peut 
se faire soit au moyen d’une approche purement numérique, soit au moyen de techniques 
semi-analytique. La réalisation de ce projet a comporté les deux étapes clés suivantes: 

En premier, nous avons utilisé la méthode de la transformée de Fourier à pas divisé 
pour simuler et analyser l’influence de l’effet de perte sur la dynamique de propagation des 
solitons optiques. Nous nous sommes retrouvés devant une amplitude du soliton qui diminue 
énormément le long de la fibre en présence de perte et en son absence, le soliton reste stable 
et se propage sans distorsion. Les pertes étant un facteur limitatif qui réduit la puissance du 
soliton et par conséquent limite les performances de la transmission optique par soliton sur les 
longues distances et ceci confirme la nécessité de l’amplification optique localisée afin de 
compenser l’atténuation de la fibre dans les télécommunications optique. Ensuite, nous nous 
sommes intéressés à l’effet de l’interaction entre les solitons fondamentaux et aussi à 
l’influence de différents paramètres. Clairement présentée au chapitre 3, l’interaction des 
solitons dépend fortement des deux paramètres qui sont l’amplitude relative et la phase. 
D’après les figures, deux solitons proches pouvaient interagir fortement d’où on assiste au 
phénomène de déformation par attirance ou répulsion. Pour limiter cette interaction qui est 
néfaste pour les transmissions optique, il est nécessaire d’espacer fortement les impulsions 
successives, ce qui conduit à une réduction du débit. Nous remarquons que l’effet d’une 
différence d’amplitude entre les deux solitons en phase présentait une propagation en 
parallèle sans collision ni répulsion des deux solitons. 
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L’évolution des solitons d’ordre supérieurs a été aussi étudiée et elle a montré clairement une 
succession de collision et d’éclatements. 

Dans la deuxième partie, nous avons étudié analytiquement l’existence des ondes 
solitaires dans un réseau de Bragg fibré en régime femtosecondes, où la propagation des 
impulsions est décrite par un système d’équations des modes couplés incluant des termes 
d’ordre élevé, des non linéarités de troisième, cinquième et septième ordre Afin de résoudre 
ce système d’équations nous avons utilisé en premier la méthode perturbative multi-échelle et 
le système d’équations non linéaires couplées (NLCME) s’est réduit en une seule équation 
dite l’équation de Schrôdinger non linéaire perturbée (PNLS) qui décrit la propagation non 
linéaire des impulsions aux bords de la zone du gap de la bande photonique. Ensuite, nous 
avons appliqué la fonction de l’amplitude complexe combinée pour résoudre cette équation 
compliqué qui contient des termes non-Kerr et des termes de dispersion de troisième et 
quatrième ordre. Les solutions obtenues sont de type solitons noirs et solitons brillants. En 
outre, nous avons également présenté les conditions d’existences de ces solutions et leurs 
paramètres physiques importants tels que la puissance et la largeur d’impulsion. 

Finalement, les travaux de recherche entrepris dans le cadre de cette thèse, peuvent envisager 
de nombreuses perspectives, qui méritent à notre avis d’être approfondies et ceci concerne les 
aspects théoriques suivants : 

• L’étude de la stabilité des solitons Bragg, 

• L’étude de l’interaction entre les solitons Bragg 

• Le développement d’une technique numérique pour résoudre les équations des modes 
couplés non-linéaires à ordre élevé. 
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